library(ggplot2)
library(dplyr)
## 
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
library(tidyverse)
## ── Attaching packages ──────────────────────────────────────────────── tidyverse 1.3.0 ──
## ✓ tibble  3.0.0     ✓ purrr   0.3.4
## ✓ tidyr   1.0.2     ✓ stringr 1.4.0
## ✓ readr   1.3.1     ✓ forcats 0.5.0
## ── Conflicts ─────────────────────────────────────────────────── tidyverse_conflicts() ──
## x dplyr::filter() masks stats::filter()
## x dplyr::lag()    masks stats::lag()
library(knitr)
library(gridExtra)
## 
## Attaching package: 'gridExtra'
## The following object is masked from 'package:dplyr':
## 
##     combine
library(stringr)
library(purrr)
samples_sinus <- read_csv("samples_normal_sinus_rythmn.csv", skip = 1)
## Parsed with column specification:
## cols(
##   `'seconds'` = col_double(),
##   `'mV'` = col_double()
## )
colnames(samples_sinus) <- c("time", "signal")
S <- 5000
burn <- 5000
y <- samples_sinus$signal[55:210]
t <- samples_sinus$time[55:210]
beta_curr <- 1
alpha_curr <- rep(1, 12)
delta_curr<- seq(0, max(t), length.out = 17) 
end_time <- max(t)
tao_curr <- 1
mean_curr <- rep(mean(y), length(y)) #mean_curr is a vector cause mean_star of y is a vector so I think r only lets me run it if it is a vector
dc_curr <- 0.5
miu_curr <- mean_curr - dc_curr

#storage

DC <- NULL
ALPHA <- rep(list(NULL), 12)
BETA <- NULL
DELTA <- rep(list(NULL), 17)
TAO <- NULL
MEAN <- rep(list(NULL), length(y))
for(s in 1:S) {
  print(c("s: ", s))
  
  #1 draw of Beta
  
  beta_lst <- beta_one_samp(beta_curr, BETA, s)
  beta_curr <- beta_lst[[1]]
  BETA <- beta_lst[[2]]
  #1 draw of Alphas (1-12)
  
  for (k in 1:12) {
    alpha_lst <- alpha_one_samp(k, alpha_curr[k], ALPHA[[k]], s)
    alpha_curr[k] <- alpha_lst[[1]]
    ALPHA[[k]] <- alpha_lst[[2]]
  }
  
  #1 draw of Deltas (1-17)
  
  delta_lst <- delta_one_samp(delta_curr[1], 0, delta_curr[2], DELTA[[1]], s)
  delta_curr[1] <- delta_lst[[1]]
  DELTA[[1]] <- delta_lst[[2]]
  
  for (k in 2:16) {
    delta_lst <- delta_one_samp(delta_curr[k], delta_curr[k-1], delta_curr[k+1], DELTA[[k]], s)
    delta_curr[k] <- delta_lst[[1]]
    DELTA[[k]] <- delta_lst[[2]]
  }
  
  delta_lst <- delta_one_samp(delta_curr[17], delta_curr[[16]], end_time,  DELTA[[17]], s)
  delta_curr[17] <- delta_lst[[1]]
  DELTA[[17]] <- delta_lst[[2]]
  
  #1 draw of Tao 
  tao_lst <- tao_one_samp_cond(tao_curr, TAO, s)
  tao_curr <- tao_lst[[1]]
  TAO <- tao_lst[[2]]

  #1 draw of DC 
  dc_lst <- dc_one_samp(dc_curr, DC, s)
  dc_curr <- dc_lst[[1]]
  DC <- dc_lst[[2]]
  
  mean_curr <- dc_curr + gen_miu_star(t, alpha_curr = alpha_curr, beta_curr = beta_curr, delta_curr = delta_curr)
  for (k in 1:length(y)) {
    MEAN[[k]] <- c(MEAN[[k]], mean_curr[[k]])
  }
  
  #print(c("MEAN CURR: ", mean_curr))
}
## [1] "s: " "1"  
## [1] "s: " "2"  
## [1] "s: " "3"  
## [1] "s: " "4"  
## [1] "s: " "5"  
## [1] "s: " "6"  
## [1] "s: " "7"  
## [1] "s: " "8"  
## [1] "s: " "9"  
## [1] "s: " "10" 
## [1] "s: " "11" 
## [1] "s: " "12" 
## [1] "s: " "13" 
## [1] "s: " "14" 
## [1] "s: " "15" 
## [1] "s: " "16" 
## [1] "s: " "17" 
## [1] "s: " "18" 
## [1] "s: " "19" 
## [1] "s: " "20" 
## [1] "s: " "21" 
## [1] "s: " "22" 
## [1] "s: " "23" 
## [1] "s: " "24" 
## [1] "s: " "25" 
## [1] "s: " "26" 
## [1] "s: " "27" 
## [1] "s: " "28" 
## [1] "s: " "29" 
## [1] "s: " "30" 
## [1] "s: " "31" 
## [1] "s: " "32" 
## [1] "s: " "33" 
## [1] "s: " "34" 
## [1] "s: " "35" 
## [1] "s: " "36" 
## [1] "s: " "37" 
## [1] "s: " "38" 
## [1] "s: " "39" 
## [1] "s: " "40" 
## [1] "s: " "41" 
## [1] "s: " "42" 
## [1] "s: " "43" 
## [1] "s: " "44" 
## [1] "s: " "45" 
## [1] "s: " "46" 
## [1] "s: " "47" 
## [1] "s: " "48" 
## [1] "s: " "49" 
## [1] "s: " "50" 
## [1] "s: " "51" 
## [1] "s: " "52" 
## [1] "s: " "53" 
## [1] "s: " "54" 
## [1] "s: " "55" 
## [1] "s: " "56" 
## [1] "s: " "57" 
## [1] "s: " "58" 
## [1] "s: " "59" 
## [1] "s: " "60" 
## [1] "s: " "61" 
## [1] "s: " "62" 
## [1] "s: " "63" 
## [1] "s: " "64" 
## [1] "s: " "65" 
## [1] "s: " "66" 
## [1] "s: " "67" 
## [1] "s: " "68" 
## [1] "s: " "69" 
## [1] "s: " "70" 
## [1] "s: " "71" 
## [1] "s: " "72" 
## [1] "s: " "73" 
## [1] "s: " "74" 
## [1] "s: " "75" 
## [1] "s: " "76" 
## [1] "s: " "77" 
## [1] "s: " "78" 
## [1] "s: " "79" 
## [1] "s: " "80" 
## [1] "s: " "81" 
## [1] "s: " "82" 
## [1] "s: " "83" 
## [1] "s: " "84" 
## [1] "s: " "85" 
## [1] "s: " "86" 
## [1] "s: " "87" 
## [1] "s: " "88" 
## [1] "s: " "89" 
## [1] "s: " "90" 
## [1] "s: " "91" 
## [1] "s: " "92" 
## [1] "s: " "93" 
## [1] "s: " "94" 
## [1] "s: " "95" 
## [1] "s: " "96" 
## [1] "s: " "97" 
## [1] "s: " "98" 
## [1] "s: " "99" 
## [1] "s: " "100"
## [1] "s: " "101"
## [1] "s: " "102"
## [1] "s: " "103"
## [1] "s: " "104"
## [1] "s: " "105"
## [1] "s: " "106"
## [1] "s: " "107"
## [1] "s: " "108"
## [1] "s: " "109"
## [1] "s: " "110"
## [1] "s: " "111"
## [1] "s: " "112"
## [1] "s: " "113"
## [1] "s: " "114"
## [1] "s: " "115"
## [1] "s: " "116"
## [1] "s: " "117"
## [1] "s: " "118"
## [1] "s: " "119"
## [1] "s: " "120"
## [1] "s: " "121"
## [1] "s: " "122"
## [1] "s: " "123"
## [1] "s: " "124"
## [1] "s: " "125"
## [1] "s: " "126"
## [1] "s: " "127"
## [1] "s: " "128"
## [1] "s: " "129"
## [1] "s: " "130"
## [1] "s: " "131"
## [1] "s: " "132"
## [1] "s: " "133"
## [1] "s: " "134"
## [1] "s: " "135"
## [1] "s: " "136"
## [1] "s: " "137"
## [1] "s: " "138"
## [1] "s: " "139"
## [1] "s: " "140"
## [1] "s: " "141"
## [1] "s: " "142"
## [1] "s: " "143"
## [1] "s: " "144"
## [1] "s: " "145"
## [1] "s: " "146"
## [1] "s: " "147"
## [1] "s: " "148"
## [1] "s: " "149"
## [1] "s: " "150"
## [1] "s: " "151"
## [1] "s: " "152"
## [1] "s: " "153"
## [1] "s: " "154"
## [1] "s: " "155"
## [1] "s: " "156"
## [1] "s: " "157"
## [1] "s: " "158"
## [1] "s: " "159"
## [1] "s: " "160"
## [1] "s: " "161"
## [1] "s: " "162"
## [1] "s: " "163"
## [1] "s: " "164"
## [1] "s: " "165"
## [1] "s: " "166"
## [1] "s: " "167"
## [1] "s: " "168"
## [1] "s: " "169"
## [1] "s: " "170"
## [1] "s: " "171"
## [1] "s: " "172"
## [1] "s: " "173"
## [1] "s: " "174"
## [1] "s: " "175"
## [1] "s: " "176"
## [1] "s: " "177"
## [1] "s: " "178"
## [1] "s: " "179"
## [1] "s: " "180"
## [1] "s: " "181"
## [1] "s: " "182"
## [1] "s: " "183"
## [1] "s: " "184"
## [1] "s: " "185"
## [1] "s: " "186"
## [1] "s: " "187"
## [1] "s: " "188"
## [1] "s: " "189"
## [1] "s: " "190"
## [1] "s: " "191"
## [1] "s: " "192"
## [1] "s: " "193"
## [1] "s: " "194"
## [1] "s: " "195"
## [1] "s: " "196"
## [1] "s: " "197"
## [1] "s: " "198"
## [1] "s: " "199"
## [1] "s: " "200"
## [1] "s: " "201"
## [1] "s: " "202"
## [1] "s: " "203"
## [1] "s: " "204"
## [1] "s: " "205"
## [1] "s: " "206"
## [1] "s: " "207"
## [1] "s: " "208"
## [1] "s: " "209"
## [1] "s: " "210"
## [1] "s: " "211"
## [1] "s: " "212"
## [1] "s: " "213"
## [1] "s: " "214"
## [1] "s: " "215"
## [1] "s: " "216"
## [1] "s: " "217"
## [1] "s: " "218"
## [1] "s: " "219"
## [1] "s: " "220"
## [1] "s: " "221"
## [1] "s: " "222"
## [1] "s: " "223"
## [1] "s: " "224"
## [1] "s: " "225"
## [1] "s: " "226"
## [1] "s: " "227"
## [1] "s: " "228"
## [1] "s: " "229"
## [1] "s: " "230"
## [1] "s: " "231"
## [1] "s: " "232"
## [1] "s: " "233"
## [1] "s: " "234"
## [1] "s: " "235"
## [1] "s: " "236"
## [1] "s: " "237"
## [1] "s: " "238"
## [1] "s: " "239"
## [1] "s: " "240"
## [1] "s: " "241"
## [1] "s: " "242"
## [1] "s: " "243"
## [1] "s: " "244"
## [1] "s: " "245"
## [1] "s: " "246"
## [1] "s: " "247"
## [1] "s: " "248"
## [1] "s: " "249"
## [1] "s: " "250"
## [1] "s: " "251"
## [1] "s: " "252"
## [1] "s: " "253"
## [1] "s: " "254"
## [1] "s: " "255"
## [1] "s: " "256"
## [1] "s: " "257"
## [1] "s: " "258"
## [1] "s: " "259"
## [1] "s: " "260"
## [1] "s: " "261"
## [1] "s: " "262"
## [1] "s: " "263"
## [1] "s: " "264"
## [1] "s: " "265"
## [1] "s: " "266"
## [1] "s: " "267"
## [1] "s: " "268"
## [1] "s: " "269"
## [1] "s: " "270"
## [1] "s: " "271"
## [1] "s: " "272"
## [1] "s: " "273"
## [1] "s: " "274"
## [1] "s: " "275"
## [1] "s: " "276"
## [1] "s: " "277"
## [1] "s: " "278"
## [1] "s: " "279"
## [1] "s: " "280"
## [1] "s: " "281"
## [1] "s: " "282"
## [1] "s: " "283"
## [1] "s: " "284"
## [1] "s: " "285"
## [1] "s: " "286"
## [1] "s: " "287"
## [1] "s: " "288"
## [1] "s: " "289"
## [1] "s: " "290"
## [1] "s: " "291"
## [1] "s: " "292"
## [1] "s: " "293"
## [1] "s: " "294"
## [1] "s: " "295"
## [1] "s: " "296"
## [1] "s: " "297"
## [1] "s: " "298"
## [1] "s: " "299"
## [1] "s: " "300"
## [1] "s: " "301"
## [1] "s: " "302"
## [1] "s: " "303"
## [1] "s: " "304"
## [1] "s: " "305"
## [1] "s: " "306"
## [1] "s: " "307"
## [1] "s: " "308"
## [1] "s: " "309"
## [1] "s: " "310"
## [1] "s: " "311"
## [1] "s: " "312"
## [1] "s: " "313"
## [1] "s: " "314"
## [1] "s: " "315"
## [1] "s: " "316"
## [1] "s: " "317"
## [1] "s: " "318"
## [1] "s: " "319"
## [1] "s: " "320"
## [1] "s: " "321"
## [1] "s: " "322"
## [1] "s: " "323"
## [1] "s: " "324"
## [1] "s: " "325"
## [1] "s: " "326"
## [1] "s: " "327"
## [1] "s: " "328"
## [1] "s: " "329"
## [1] "s: " "330"
## [1] "s: " "331"
## [1] "s: " "332"
## [1] "s: " "333"
## [1] "s: " "334"
## [1] "s: " "335"
## [1] "s: " "336"
## [1] "s: " "337"
## [1] "s: " "338"
## [1] "s: " "339"
## [1] "s: " "340"
## [1] "s: " "341"
## [1] "s: " "342"
## [1] "s: " "343"
## [1] "s: " "344"
## [1] "s: " "345"
## [1] "s: " "346"
## [1] "s: " "347"
## [1] "s: " "348"
## [1] "s: " "349"
## [1] "s: " "350"
## [1] "s: " "351"
## [1] "s: " "352"
## [1] "s: " "353"
## [1] "s: " "354"
## [1] "s: " "355"
## [1] "s: " "356"
## [1] "s: " "357"
## [1] "s: " "358"
## [1] "s: " "359"
## [1] "s: " "360"
## [1] "s: " "361"
## [1] "s: " "362"
## [1] "s: " "363"
## [1] "s: " "364"
## [1] "s: " "365"
## [1] "s: " "366"
## [1] "s: " "367"
## [1] "s: " "368"
## [1] "s: " "369"
## [1] "s: " "370"
## [1] "s: " "371"
## [1] "s: " "372"
## [1] "s: " "373"
## [1] "s: " "374"
## [1] "s: " "375"
## [1] "s: " "376"
## [1] "s: " "377"
## [1] "s: " "378"
## [1] "s: " "379"
## [1] "s: " "380"
## [1] "s: " "381"
## [1] "s: " "382"
## [1] "s: " "383"
## [1] "s: " "384"
## [1] "s: " "385"
## [1] "s: " "386"
## [1] "s: " "387"
## [1] "s: " "388"
## [1] "s: " "389"
## [1] "s: " "390"
## [1] "s: " "391"
## [1] "s: " "392"
## [1] "s: " "393"
## [1] "s: " "394"
## [1] "s: " "395"
## [1] "s: " "396"
## [1] "s: " "397"
## [1] "s: " "398"
## [1] "s: " "399"
## [1] "s: " "400"
## [1] "s: " "401"
## [1] "s: " "402"
## [1] "s: " "403"
## [1] "s: " "404"
## [1] "s: " "405"
## [1] "s: " "406"
## [1] "s: " "407"
## [1] "s: " "408"
## [1] "s: " "409"
## [1] "s: " "410"
## [1] "s: " "411"
## [1] "s: " "412"
## [1] "s: " "413"
## [1] "s: " "414"
## [1] "s: " "415"
## [1] "s: " "416"
## [1] "s: " "417"
## [1] "s: " "418"
## [1] "s: " "419"
## [1] "s: " "420"
## [1] "s: " "421"
## [1] "s: " "422"
## [1] "s: " "423"
## [1] "s: " "424"
## [1] "s: " "425"
## [1] "s: " "426"
## [1] "s: " "427"
## [1] "s: " "428"
## [1] "s: " "429"
## [1] "s: " "430"
## [1] "s: " "431"
## [1] "s: " "432"
## [1] "s: " "433"
## [1] "s: " "434"
## [1] "s: " "435"
## [1] "s: " "436"
## [1] "s: " "437"
## [1] "s: " "438"
## [1] "s: " "439"
## [1] "s: " "440"
## [1] "s: " "441"
## [1] "s: " "442"
## [1] "s: " "443"
## [1] "s: " "444"
## [1] "s: " "445"
## [1] "s: " "446"
## [1] "s: " "447"
## [1] "s: " "448"
## [1] "s: " "449"
## [1] "s: " "450"
## [1] "s: " "451"
## [1] "s: " "452"
## [1] "s: " "453"
## [1] "s: " "454"
## [1] "s: " "455"
## [1] "s: " "456"
## [1] "s: " "457"
## [1] "s: " "458"
## [1] "s: " "459"
## [1] "s: " "460"
## [1] "s: " "461"
## [1] "s: " "462"
## [1] "s: " "463"
## [1] "s: " "464"
## [1] "s: " "465"
## [1] "s: " "466"
## [1] "s: " "467"
## [1] "s: " "468"
## [1] "s: " "469"
## [1] "s: " "470"
## [1] "s: " "471"
## [1] "s: " "472"
## [1] "s: " "473"
## [1] "s: " "474"
## [1] "s: " "475"
## [1] "s: " "476"
## [1] "s: " "477"
## [1] "s: " "478"
## [1] "s: " "479"
## [1] "s: " "480"
## [1] "s: " "481"
## [1] "s: " "482"
## [1] "s: " "483"
## [1] "s: " "484"
## [1] "s: " "485"
## [1] "s: " "486"
## [1] "s: " "487"
## [1] "s: " "488"
## [1] "s: " "489"
## [1] "s: " "490"
## [1] "s: " "491"
## [1] "s: " "492"
## [1] "s: " "493"
## [1] "s: " "494"
## [1] "s: " "495"
## [1] "s: " "496"
## [1] "s: " "497"
## [1] "s: " "498"
## [1] "s: " "499"
## [1] "s: " "500"
## [1] "s: " "501"
## [1] "s: " "502"
## [1] "s: " "503"
## [1] "s: " "504"
## [1] "s: " "505"
## [1] "s: " "506"
## [1] "s: " "507"
## [1] "s: " "508"
## [1] "s: " "509"
## [1] "s: " "510"
## [1] "s: " "511"
## [1] "s: " "512"
## [1] "s: " "513"
## [1] "s: " "514"
## [1] "s: " "515"
## [1] "s: " "516"
## [1] "s: " "517"
## [1] "s: " "518"
## [1] "s: " "519"
## [1] "s: " "520"
## [1] "s: " "521"
## [1] "s: " "522"
## [1] "s: " "523"
## [1] "s: " "524"
## [1] "s: " "525"
## [1] "s: " "526"
## [1] "s: " "527"
## [1] "s: " "528"
## [1] "s: " "529"
## [1] "s: " "530"
## [1] "s: " "531"
## [1] "s: " "532"
## [1] "s: " "533"
## [1] "s: " "534"
## [1] "s: " "535"
## [1] "s: " "536"
## [1] "s: " "537"
## [1] "s: " "538"
## [1] "s: " "539"
## [1] "s: " "540"
## [1] "s: " "541"
## [1] "s: " "542"
## [1] "s: " "543"
## [1] "s: " "544"
## [1] "s: " "545"
## [1] "s: " "546"
## [1] "s: " "547"
## [1] "s: " "548"
## [1] "s: " "549"
## [1] "s: " "550"
## [1] "s: " "551"
## [1] "s: " "552"
## [1] "s: " "553"
## [1] "s: " "554"
## [1] "s: " "555"
## [1] "s: " "556"
## [1] "s: " "557"
## [1] "s: " "558"
## [1] "s: " "559"
## [1] "s: " "560"
## [1] "s: " "561"
## [1] "s: " "562"
## [1] "s: " "563"
## [1] "s: " "564"
## [1] "s: " "565"
## [1] "s: " "566"
## [1] "s: " "567"
## [1] "s: " "568"
## [1] "s: " "569"
## [1] "s: " "570"
## [1] "s: " "571"
## [1] "s: " "572"
## [1] "s: " "573"
## [1] "s: " "574"
## [1] "s: " "575"
## [1] "s: " "576"
## [1] "s: " "577"
## [1] "s: " "578"
## [1] "s: " "579"
## [1] "s: " "580"
## [1] "s: " "581"
## [1] "s: " "582"
## [1] "s: " "583"
## [1] "s: " "584"
## [1] "s: " "585"
## [1] "s: " "586"
## [1] "s: " "587"
## [1] "s: " "588"
## [1] "s: " "589"
## [1] "s: " "590"
## [1] "s: " "591"
## [1] "s: " "592"
## [1] "s: " "593"
## [1] "s: " "594"
## [1] "s: " "595"
## [1] "s: " "596"
## [1] "s: " "597"
## [1] "s: " "598"
## [1] "s: " "599"
## [1] "s: " "600"
## [1] "s: " "601"
## [1] "s: " "602"
## [1] "s: " "603"
## [1] "s: " "604"
## [1] "s: " "605"
## [1] "s: " "606"
## [1] "s: " "607"
## [1] "s: " "608"
## [1] "s: " "609"
## [1] "s: " "610"
## [1] "s: " "611"
## [1] "s: " "612"
## [1] "s: " "613"
## [1] "s: " "614"
## [1] "s: " "615"
## [1] "s: " "616"
## [1] "s: " "617"
## [1] "s: " "618"
## [1] "s: " "619"
## [1] "s: " "620"
## [1] "s: " "621"
## [1] "s: " "622"
## [1] "s: " "623"
## [1] "s: " "624"
## [1] "s: " "625"
## [1] "s: " "626"
## [1] "s: " "627"
## [1] "s: " "628"
## [1] "s: " "629"
## [1] "s: " "630"
## [1] "s: " "631"
## [1] "s: " "632"
## [1] "s: " "633"
## [1] "s: " "634"
## [1] "s: " "635"
## [1] "s: " "636"
## [1] "s: " "637"
## [1] "s: " "638"
## [1] "s: " "639"
## [1] "s: " "640"
## [1] "s: " "641"
## [1] "s: " "642"
## [1] "s: " "643"
## [1] "s: " "644"
## [1] "s: " "645"
## [1] "s: " "646"
## [1] "s: " "647"
## [1] "s: " "648"
## [1] "s: " "649"
## [1] "s: " "650"
## [1] "s: " "651"
## [1] "s: " "652"
## [1] "s: " "653"
## [1] "s: " "654"
## [1] "s: " "655"
## [1] "s: " "656"
## [1] "s: " "657"
## [1] "s: " "658"
## [1] "s: " "659"
## [1] "s: " "660"
## [1] "s: " "661"
## [1] "s: " "662"
## [1] "s: " "663"
## [1] "s: " "664"
## [1] "s: " "665"
## [1] "s: " "666"
## [1] "s: " "667"
## [1] "s: " "668"
## [1] "s: " "669"
## [1] "s: " "670"
## [1] "s: " "671"
## [1] "s: " "672"
## [1] "s: " "673"
## [1] "s: " "674"
## [1] "s: " "675"
## [1] "s: " "676"
## [1] "s: " "677"
## [1] "s: " "678"
## [1] "s: " "679"
## [1] "s: " "680"
## [1] "s: " "681"
## [1] "s: " "682"
## [1] "s: " "683"
## [1] "s: " "684"
## [1] "s: " "685"
## [1] "s: " "686"
## [1] "s: " "687"
## [1] "s: " "688"
## [1] "s: " "689"
## [1] "s: " "690"
## [1] "s: " "691"
## [1] "s: " "692"
## [1] "s: " "693"
## [1] "s: " "694"
## [1] "s: " "695"
## [1] "s: " "696"
## [1] "s: " "697"
## [1] "s: " "698"
## [1] "s: " "699"
## [1] "s: " "700"
## [1] "s: " "701"
## [1] "s: " "702"
## [1] "s: " "703"
## [1] "s: " "704"
## [1] "s: " "705"
## [1] "s: " "706"
## [1] "s: " "707"
## [1] "s: " "708"
## [1] "s: " "709"
## [1] "s: " "710"
## [1] "s: " "711"
## [1] "s: " "712"
## [1] "s: " "713"
## [1] "s: " "714"
## [1] "s: " "715"
## [1] "s: " "716"
## [1] "s: " "717"
## [1] "s: " "718"
## [1] "s: " "719"
## [1] "s: " "720"
## [1] "s: " "721"
## [1] "s: " "722"
## [1] "s: " "723"
## [1] "s: " "724"
## [1] "s: " "725"
## [1] "s: " "726"
## [1] "s: " "727"
## [1] "s: " "728"
## [1] "s: " "729"
## [1] "s: " "730"
## [1] "s: " "731"
## [1] "s: " "732"
## [1] "s: " "733"
## [1] "s: " "734"
## [1] "s: " "735"
## [1] "s: " "736"
## [1] "s: " "737"
## [1] "s: " "738"
## [1] "s: " "739"
## [1] "s: " "740"
## [1] "s: " "741"
## [1] "s: " "742"
## [1] "s: " "743"
## [1] "s: " "744"
## [1] "s: " "745"
## [1] "s: " "746"
## [1] "s: " "747"
## [1] "s: " "748"
## [1] "s: " "749"
## [1] "s: " "750"
## [1] "s: " "751"
## [1] "s: " "752"
## [1] "s: " "753"
## [1] "s: " "754"
## [1] "s: " "755"
## [1] "s: " "756"
## [1] "s: " "757"
## [1] "s: " "758"
## [1] "s: " "759"
## [1] "s: " "760"
## [1] "s: " "761"
## [1] "s: " "762"
## [1] "s: " "763"
## [1] "s: " "764"
## [1] "s: " "765"
## [1] "s: " "766"
## [1] "s: " "767"
## [1] "s: " "768"
## [1] "s: " "769"
## [1] "s: " "770"
## [1] "s: " "771"
## [1] "s: " "772"
## [1] "s: " "773"
## [1] "s: " "774"
## [1] "s: " "775"
## [1] "s: " "776"
## [1] "s: " "777"
## [1] "s: " "778"
## [1] "s: " "779"
## [1] "s: " "780"
## [1] "s: " "781"
## [1] "s: " "782"
## [1] "s: " "783"
## [1] "s: " "784"
## [1] "s: " "785"
## [1] "s: " "786"
## [1] "s: " "787"
## [1] "s: " "788"
## [1] "s: " "789"
## [1] "s: " "790"
## [1] "s: " "791"
## [1] "s: " "792"
## [1] "s: " "793"
## [1] "s: " "794"
## [1] "s: " "795"
## [1] "s: " "796"
## [1] "s: " "797"
## [1] "s: " "798"
## [1] "s: " "799"
## [1] "s: " "800"
## [1] "s: " "801"
## [1] "s: " "802"
## [1] "s: " "803"
## [1] "s: " "804"
## [1] "s: " "805"
## [1] "s: " "806"
## [1] "s: " "807"
## [1] "s: " "808"
## [1] "s: " "809"
## [1] "s: " "810"
## [1] "s: " "811"
## [1] "s: " "812"
## [1] "s: " "813"
## [1] "s: " "814"
## [1] "s: " "815"
## [1] "s: " "816"
## [1] "s: " "817"
## [1] "s: " "818"
## [1] "s: " "819"
## [1] "s: " "820"
## [1] "s: " "821"
## [1] "s: " "822"
## [1] "s: " "823"
## [1] "s: " "824"
## [1] "s: " "825"
## [1] "s: " "826"
## [1] "s: " "827"
## [1] "s: " "828"
## [1] "s: " "829"
## [1] "s: " "830"
## [1] "s: " "831"
## [1] "s: " "832"
## [1] "s: " "833"
## [1] "s: " "834"
## [1] "s: " "835"
## [1] "s: " "836"
## [1] "s: " "837"
## [1] "s: " "838"
## [1] "s: " "839"
## [1] "s: " "840"
## [1] "s: " "841"
## [1] "s: " "842"
## [1] "s: " "843"
## [1] "s: " "844"
## [1] "s: " "845"
## [1] "s: " "846"
## [1] "s: " "847"
## [1] "s: " "848"
## [1] "s: " "849"
## [1] "s: " "850"
## [1] "s: " "851"
## [1] "s: " "852"
## [1] "s: " "853"
## [1] "s: " "854"
## [1] "s: " "855"
## [1] "s: " "856"
## [1] "s: " "857"
## [1] "s: " "858"
## [1] "s: " "859"
## [1] "s: " "860"
## [1] "s: " "861"
## [1] "s: " "862"
## [1] "s: " "863"
## [1] "s: " "864"
## [1] "s: " "865"
## [1] "s: " "866"
## [1] "s: " "867"
## [1] "s: " "868"
## [1] "s: " "869"
## [1] "s: " "870"
## [1] "s: " "871"
## [1] "s: " "872"
## [1] "s: " "873"
## [1] "s: " "874"
## [1] "s: " "875"
## [1] "s: " "876"
## [1] "s: " "877"
## [1] "s: " "878"
## [1] "s: " "879"
## [1] "s: " "880"
## [1] "s: " "881"
## [1] "s: " "882"
## [1] "s: " "883"
## [1] "s: " "884"
## [1] "s: " "885"
## [1] "s: " "886"
## [1] "s: " "887"
## [1] "s: " "888"
## [1] "s: " "889"
## [1] "s: " "890"
## [1] "s: " "891"
## [1] "s: " "892"
## [1] "s: " "893"
## [1] "s: " "894"
## [1] "s: " "895"
## [1] "s: " "896"
## [1] "s: " "897"
## [1] "s: " "898"
## [1] "s: " "899"
## [1] "s: " "900"
## [1] "s: " "901"
## [1] "s: " "902"
## [1] "s: " "903"
## [1] "s: " "904"
## [1] "s: " "905"
## [1] "s: " "906"
## [1] "s: " "907"
## [1] "s: " "908"
## [1] "s: " "909"
## [1] "s: " "910"
## [1] "s: " "911"
## [1] "s: " "912"
## [1] "s: " "913"
## [1] "s: " "914"
## [1] "s: " "915"
## [1] "s: " "916"
## [1] "s: " "917"
## [1] "s: " "918"
## [1] "s: " "919"
## [1] "s: " "920"
## [1] "s: " "921"
## [1] "s: " "922"
## [1] "s: " "923"
## [1] "s: " "924"
## [1] "s: " "925"
## [1] "s: " "926"
## [1] "s: " "927"
## [1] "s: " "928"
## [1] "s: " "929"
## [1] "s: " "930"
## [1] "s: " "931"
## [1] "s: " "932"
## [1] "s: " "933"
## [1] "s: " "934"
## [1] "s: " "935"
## [1] "s: " "936"
## [1] "s: " "937"
## [1] "s: " "938"
## [1] "s: " "939"
## [1] "s: " "940"
## [1] "s: " "941"
## [1] "s: " "942"
## [1] "s: " "943"
## [1] "s: " "944"
## [1] "s: " "945"
## [1] "s: " "946"
## [1] "s: " "947"
## [1] "s: " "948"
## [1] "s: " "949"
## [1] "s: " "950"
## [1] "s: " "951"
## [1] "s: " "952"
## [1] "s: " "953"
## [1] "s: " "954"
## [1] "s: " "955"
## [1] "s: " "956"
## [1] "s: " "957"
## [1] "s: " "958"
## [1] "s: " "959"
## [1] "s: " "960"
## [1] "s: " "961"
## [1] "s: " "962"
## [1] "s: " "963"
## [1] "s: " "964"
## [1] "s: " "965"
## [1] "s: " "966"
## [1] "s: " "967"
## [1] "s: " "968"
## [1] "s: " "969"
## [1] "s: " "970"
## [1] "s: " "971"
## [1] "s: " "972"
## [1] "s: " "973"
## [1] "s: " "974"
## [1] "s: " "975"
## [1] "s: " "976"
## [1] "s: " "977"
## [1] "s: " "978"
## [1] "s: " "979"
## [1] "s: " "980"
## [1] "s: " "981"
## [1] "s: " "982"
## [1] "s: " "983"
## [1] "s: " "984"
## [1] "s: " "985"
## [1] "s: " "986"
## [1] "s: " "987"
## [1] "s: " "988"
## [1] "s: " "989"
## [1] "s: " "990"
## [1] "s: " "991"
## [1] "s: " "992"
## [1] "s: " "993"
## [1] "s: " "994"
## [1] "s: " "995"
## [1] "s: " "996"
## [1] "s: " "997"
## [1] "s: " "998"
## [1] "s: " "999"
## [1] "s: "  "1000"
## [1] "s: "  "1001"
## [1] "s: "  "1002"
## [1] "s: "  "1003"
## [1] "s: "  "1004"
## [1] "s: "  "1005"
## [1] "s: "  "1006"
## [1] "s: "  "1007"
## [1] "s: "  "1008"
## [1] "s: "  "1009"
## [1] "s: "  "1010"
## [1] "s: "  "1011"
## [1] "s: "  "1012"
## [1] "s: "  "1013"
## [1] "s: "  "1014"
## [1] "s: "  "1015"
## [1] "s: "  "1016"
## [1] "s: "  "1017"
## [1] "s: "  "1018"
## [1] "s: "  "1019"
## [1] "s: "  "1020"
## [1] "s: "  "1021"
## [1] "s: "  "1022"
## [1] "s: "  "1023"
## [1] "s: "  "1024"
## [1] "s: "  "1025"
## [1] "s: "  "1026"
## [1] "s: "  "1027"
## [1] "s: "  "1028"
## [1] "s: "  "1029"
## [1] "s: "  "1030"
## [1] "s: "  "1031"
## [1] "s: "  "1032"
## [1] "s: "  "1033"
## [1] "s: "  "1034"
## [1] "s: "  "1035"
## [1] "s: "  "1036"
## [1] "s: "  "1037"
## [1] "s: "  "1038"
## [1] "s: "  "1039"
## [1] "s: "  "1040"
## [1] "s: "  "1041"
## [1] "s: "  "1042"
## [1] "s: "  "1043"
## [1] "s: "  "1044"
## [1] "s: "  "1045"
## [1] "s: "  "1046"
## [1] "s: "  "1047"
## [1] "s: "  "1048"
## [1] "s: "  "1049"
## [1] "s: "  "1050"
## [1] "s: "  "1051"
## [1] "s: "  "1052"
## [1] "s: "  "1053"
## [1] "s: "  "1054"
## [1] "s: "  "1055"
## [1] "s: "  "1056"
## [1] "s: "  "1057"
## [1] "s: "  "1058"
## [1] "s: "  "1059"
## [1] "s: "  "1060"
## [1] "s: "  "1061"
## [1] "s: "  "1062"
## [1] "s: "  "1063"
## [1] "s: "  "1064"
## [1] "s: "  "1065"
## [1] "s: "  "1066"
## [1] "s: "  "1067"
## [1] "s: "  "1068"
## [1] "s: "  "1069"
## [1] "s: "  "1070"
## [1] "s: "  "1071"
## [1] "s: "  "1072"
## [1] "s: "  "1073"
## [1] "s: "  "1074"
## [1] "s: "  "1075"
## [1] "s: "  "1076"
## [1] "s: "  "1077"
## [1] "s: "  "1078"
## [1] "s: "  "1079"
## [1] "s: "  "1080"
## [1] "s: "  "1081"
## [1] "s: "  "1082"
## [1] "s: "  "1083"
## [1] "s: "  "1084"
## [1] "s: "  "1085"
## [1] "s: "  "1086"
## [1] "s: "  "1087"
## [1] "s: "  "1088"
## [1] "s: "  "1089"
## [1] "s: "  "1090"
## [1] "s: "  "1091"
## [1] "s: "  "1092"
## [1] "s: "  "1093"
## [1] "s: "  "1094"
## [1] "s: "  "1095"
## [1] "s: "  "1096"
## [1] "s: "  "1097"
## [1] "s: "  "1098"
## [1] "s: "  "1099"
## [1] "s: "  "1100"
## [1] "s: "  "1101"
## [1] "s: "  "1102"
## [1] "s: "  "1103"
## [1] "s: "  "1104"
## [1] "s: "  "1105"
## [1] "s: "  "1106"
## [1] "s: "  "1107"
## [1] "s: "  "1108"
## [1] "s: "  "1109"
## [1] "s: "  "1110"
## [1] "s: "  "1111"
## [1] "s: "  "1112"
## [1] "s: "  "1113"
## [1] "s: "  "1114"
## [1] "s: "  "1115"
## [1] "s: "  "1116"
## [1] "s: "  "1117"
## [1] "s: "  "1118"
## [1] "s: "  "1119"
## [1] "s: "  "1120"
## [1] "s: "  "1121"
## [1] "s: "  "1122"
## [1] "s: "  "1123"
## [1] "s: "  "1124"
## [1] "s: "  "1125"
## [1] "s: "  "1126"
## [1] "s: "  "1127"
## [1] "s: "  "1128"
## [1] "s: "  "1129"
## [1] "s: "  "1130"
## [1] "s: "  "1131"
## [1] "s: "  "1132"
## [1] "s: "  "1133"
## [1] "s: "  "1134"
## [1] "s: "  "1135"
## [1] "s: "  "1136"
## [1] "s: "  "1137"
## [1] "s: "  "1138"
## [1] "s: "  "1139"
## [1] "s: "  "1140"
## [1] "s: "  "1141"
## [1] "s: "  "1142"
## [1] "s: "  "1143"
## [1] "s: "  "1144"
## [1] "s: "  "1145"
## [1] "s: "  "1146"
## [1] "s: "  "1147"
## [1] "s: "  "1148"
## [1] "s: "  "1149"
## [1] "s: "  "1150"
## [1] "s: "  "1151"
## [1] "s: "  "1152"
## [1] "s: "  "1153"
## [1] "s: "  "1154"
## [1] "s: "  "1155"
## [1] "s: "  "1156"
## [1] "s: "  "1157"
## [1] "s: "  "1158"
## [1] "s: "  "1159"
## [1] "s: "  "1160"
## [1] "s: "  "1161"
## [1] "s: "  "1162"
## [1] "s: "  "1163"
## [1] "s: "  "1164"
## [1] "s: "  "1165"
## [1] "s: "  "1166"
## [1] "s: "  "1167"
## [1] "s: "  "1168"
## [1] "s: "  "1169"
## [1] "s: "  "1170"
## [1] "s: "  "1171"
## [1] "s: "  "1172"
## [1] "s: "  "1173"
## [1] "s: "  "1174"
## [1] "s: "  "1175"
## [1] "s: "  "1176"
## [1] "s: "  "1177"
## [1] "s: "  "1178"
## [1] "s: "  "1179"
## [1] "s: "  "1180"
## [1] "s: "  "1181"
## [1] "s: "  "1182"
## [1] "s: "  "1183"
## [1] "s: "  "1184"
## [1] "s: "  "1185"
## [1] "s: "  "1186"
## [1] "s: "  "1187"
## [1] "s: "  "1188"
## [1] "s: "  "1189"
## [1] "s: "  "1190"
## [1] "s: "  "1191"
## [1] "s: "  "1192"
## [1] "s: "  "1193"
## [1] "s: "  "1194"
## [1] "s: "  "1195"
## [1] "s: "  "1196"
## [1] "s: "  "1197"
## [1] "s: "  "1198"
## [1] "s: "  "1199"
## [1] "s: "  "1200"
## [1] "s: "  "1201"
## [1] "s: "  "1202"
## [1] "s: "  "1203"
## [1] "s: "  "1204"
## [1] "s: "  "1205"
## [1] "s: "  "1206"
## [1] "s: "  "1207"
## [1] "s: "  "1208"
## [1] "s: "  "1209"
## [1] "s: "  "1210"
## [1] "s: "  "1211"
## [1] "s: "  "1212"
## [1] "s: "  "1213"
## [1] "s: "  "1214"
## [1] "s: "  "1215"
## [1] "s: "  "1216"
## [1] "s: "  "1217"
## [1] "s: "  "1218"
## [1] "s: "  "1219"
## [1] "s: "  "1220"
## [1] "s: "  "1221"
## [1] "s: "  "1222"
## [1] "s: "  "1223"
## [1] "s: "  "1224"
## [1] "s: "  "1225"
## [1] "s: "  "1226"
## [1] "s: "  "1227"
## [1] "s: "  "1228"
## [1] "s: "  "1229"
## [1] "s: "  "1230"
## [1] "s: "  "1231"
## [1] "s: "  "1232"
## [1] "s: "  "1233"
## [1] "s: "  "1234"
## [1] "s: "  "1235"
## [1] "s: "  "1236"
## [1] "s: "  "1237"
## [1] "s: "  "1238"
## [1] "s: "  "1239"
## [1] "s: "  "1240"
## [1] "s: "  "1241"
## [1] "s: "  "1242"
## [1] "s: "  "1243"
## [1] "s: "  "1244"
## [1] "s: "  "1245"
## [1] "s: "  "1246"
## [1] "s: "  "1247"
## [1] "s: "  "1248"
## [1] "s: "  "1249"
## [1] "s: "  "1250"
## [1] "s: "  "1251"
## [1] "s: "  "1252"
## [1] "s: "  "1253"
## [1] "s: "  "1254"
## [1] "s: "  "1255"
## [1] "s: "  "1256"
## [1] "s: "  "1257"
## [1] "s: "  "1258"
## [1] "s: "  "1259"
## [1] "s: "  "1260"
## [1] "s: "  "1261"
## [1] "s: "  "1262"
## [1] "s: "  "1263"
## [1] "s: "  "1264"
## [1] "s: "  "1265"
## [1] "s: "  "1266"
## [1] "s: "  "1267"
## [1] "s: "  "1268"
## [1] "s: "  "1269"
## [1] "s: "  "1270"
## [1] "s: "  "1271"
## [1] "s: "  "1272"
## [1] "s: "  "1273"
## [1] "s: "  "1274"
## [1] "s: "  "1275"
## [1] "s: "  "1276"
## [1] "s: "  "1277"
## [1] "s: "  "1278"
## [1] "s: "  "1279"
## [1] "s: "  "1280"
## [1] "s: "  "1281"
## [1] "s: "  "1282"
## [1] "s: "  "1283"
## [1] "s: "  "1284"
## [1] "s: "  "1285"
## [1] "s: "  "1286"
## [1] "s: "  "1287"
## [1] "s: "  "1288"
## [1] "s: "  "1289"
## [1] "s: "  "1290"
## [1] "s: "  "1291"
## [1] "s: "  "1292"
## [1] "s: "  "1293"
## [1] "s: "  "1294"
## [1] "s: "  "1295"
## [1] "s: "  "1296"
## [1] "s: "  "1297"
## [1] "s: "  "1298"
## [1] "s: "  "1299"
## [1] "s: "  "1300"
## [1] "s: "  "1301"
## [1] "s: "  "1302"
## [1] "s: "  "1303"
## [1] "s: "  "1304"
## [1] "s: "  "1305"
## [1] "s: "  "1306"
## [1] "s: "  "1307"
## [1] "s: "  "1308"
## [1] "s: "  "1309"
## [1] "s: "  "1310"
## [1] "s: "  "1311"
## [1] "s: "  "1312"
## [1] "s: "  "1313"
## [1] "s: "  "1314"
## [1] "s: "  "1315"
## [1] "s: "  "1316"
## [1] "s: "  "1317"
## [1] "s: "  "1318"
## [1] "s: "  "1319"
## [1] "s: "  "1320"
## [1] "s: "  "1321"
## [1] "s: "  "1322"
## [1] "s: "  "1323"
## [1] "s: "  "1324"
## [1] "s: "  "1325"
## [1] "s: "  "1326"
## [1] "s: "  "1327"
## [1] "s: "  "1328"
## [1] "s: "  "1329"
## [1] "s: "  "1330"
## [1] "s: "  "1331"
## [1] "s: "  "1332"
## [1] "s: "  "1333"
## [1] "s: "  "1334"
## [1] "s: "  "1335"
## [1] "s: "  "1336"
## [1] "s: "  "1337"
## [1] "s: "  "1338"
## [1] "s: "  "1339"
## [1] "s: "  "1340"
## [1] "s: "  "1341"
## [1] "s: "  "1342"
## [1] "s: "  "1343"
## [1] "s: "  "1344"
## [1] "s: "  "1345"
## [1] "s: "  "1346"
## [1] "s: "  "1347"
## [1] "s: "  "1348"
## [1] "s: "  "1349"
## [1] "s: "  "1350"
## [1] "s: "  "1351"
## [1] "s: "  "1352"
## [1] "s: "  "1353"
## [1] "s: "  "1354"
## [1] "s: "  "1355"
## [1] "s: "  "1356"
## [1] "s: "  "1357"
## [1] "s: "  "1358"
## [1] "s: "  "1359"
## [1] "s: "  "1360"
## [1] "s: "  "1361"
## [1] "s: "  "1362"
## [1] "s: "  "1363"
## [1] "s: "  "1364"
## [1] "s: "  "1365"
## [1] "s: "  "1366"
## [1] "s: "  "1367"
## [1] "s: "  "1368"
## [1] "s: "  "1369"
## [1] "s: "  "1370"
## [1] "s: "  "1371"
## [1] "s: "  "1372"
## [1] "s: "  "1373"
## [1] "s: "  "1374"
## [1] "s: "  "1375"
## [1] "s: "  "1376"
## [1] "s: "  "1377"
## [1] "s: "  "1378"
## [1] "s: "  "1379"
## [1] "s: "  "1380"
## [1] "s: "  "1381"
## [1] "s: "  "1382"
## [1] "s: "  "1383"
## [1] "s: "  "1384"
## [1] "s: "  "1385"
## [1] "s: "  "1386"
## [1] "s: "  "1387"
## [1] "s: "  "1388"
## [1] "s: "  "1389"
## [1] "s: "  "1390"
## [1] "s: "  "1391"
## [1] "s: "  "1392"
## [1] "s: "  "1393"
## [1] "s: "  "1394"
## [1] "s: "  "1395"
## [1] "s: "  "1396"
## [1] "s: "  "1397"
## [1] "s: "  "1398"
## [1] "s: "  "1399"
## [1] "s: "  "1400"
## [1] "s: "  "1401"
## [1] "s: "  "1402"
## [1] "s: "  "1403"
## [1] "s: "  "1404"
## [1] "s: "  "1405"
## [1] "s: "  "1406"
## [1] "s: "  "1407"
## [1] "s: "  "1408"
## [1] "s: "  "1409"
## [1] "s: "  "1410"
## [1] "s: "  "1411"
## [1] "s: "  "1412"
## [1] "s: "  "1413"
## [1] "s: "  "1414"
## [1] "s: "  "1415"
## [1] "s: "  "1416"
## [1] "s: "  "1417"
## [1] "s: "  "1418"
## [1] "s: "  "1419"
## [1] "s: "  "1420"
## [1] "s: "  "1421"
## [1] "s: "  "1422"
## [1] "s: "  "1423"
## [1] "s: "  "1424"
## [1] "s: "  "1425"
## [1] "s: "  "1426"
## [1] "s: "  "1427"
## [1] "s: "  "1428"
## [1] "s: "  "1429"
## [1] "s: "  "1430"
## [1] "s: "  "1431"
## [1] "s: "  "1432"
## [1] "s: "  "1433"
## [1] "s: "  "1434"
## [1] "s: "  "1435"
## [1] "s: "  "1436"
## [1] "s: "  "1437"
## [1] "s: "  "1438"
## [1] "s: "  "1439"
## [1] "s: "  "1440"
## [1] "s: "  "1441"
## [1] "s: "  "1442"
## [1] "s: "  "1443"
## [1] "s: "  "1444"
## [1] "s: "  "1445"
## [1] "s: "  "1446"
## [1] "s: "  "1447"
## [1] "s: "  "1448"
## [1] "s: "  "1449"
## [1] "s: "  "1450"
## [1] "s: "  "1451"
## [1] "s: "  "1452"
## [1] "s: "  "1453"
## [1] "s: "  "1454"
## [1] "s: "  "1455"
## [1] "s: "  "1456"
## [1] "s: "  "1457"
## [1] "s: "  "1458"
## [1] "s: "  "1459"
## [1] "s: "  "1460"
## [1] "s: "  "1461"
## [1] "s: "  "1462"
## [1] "s: "  "1463"
## [1] "s: "  "1464"
## [1] "s: "  "1465"
## [1] "s: "  "1466"
## [1] "s: "  "1467"
## [1] "s: "  "1468"
## [1] "s: "  "1469"
## [1] "s: "  "1470"
## [1] "s: "  "1471"
## [1] "s: "  "1472"
## [1] "s: "  "1473"
## [1] "s: "  "1474"
## [1] "s: "  "1475"
## [1] "s: "  "1476"
## [1] "s: "  "1477"
## [1] "s: "  "1478"
## [1] "s: "  "1479"
## [1] "s: "  "1480"
## [1] "s: "  "1481"
## [1] "s: "  "1482"
## [1] "s: "  "1483"
## [1] "s: "  "1484"
## [1] "s: "  "1485"
## [1] "s: "  "1486"
## [1] "s: "  "1487"
## [1] "s: "  "1488"
## [1] "s: "  "1489"
## [1] "s: "  "1490"
## [1] "s: "  "1491"
## [1] "s: "  "1492"
## [1] "s: "  "1493"
## [1] "s: "  "1494"
## [1] "s: "  "1495"
## [1] "s: "  "1496"
## [1] "s: "  "1497"
## [1] "s: "  "1498"
## [1] "s: "  "1499"
## [1] "s: "  "1500"
## [1] "s: "  "1501"
## [1] "s: "  "1502"
## [1] "s: "  "1503"
## [1] "s: "  "1504"
## [1] "s: "  "1505"
## [1] "s: "  "1506"
## [1] "s: "  "1507"
## [1] "s: "  "1508"
## [1] "s: "  "1509"
## [1] "s: "  "1510"
## [1] "s: "  "1511"
## [1] "s: "  "1512"
## [1] "s: "  "1513"
## [1] "s: "  "1514"
## [1] "s: "  "1515"
## [1] "s: "  "1516"
## [1] "s: "  "1517"
## [1] "s: "  "1518"
## [1] "s: "  "1519"
## [1] "s: "  "1520"
## [1] "s: "  "1521"
## [1] "s: "  "1522"
## [1] "s: "  "1523"
## [1] "s: "  "1524"
## [1] "s: "  "1525"
## [1] "s: "  "1526"
## [1] "s: "  "1527"
## [1] "s: "  "1528"
## [1] "s: "  "1529"
## [1] "s: "  "1530"
## [1] "s: "  "1531"
## [1] "s: "  "1532"
## [1] "s: "  "1533"
## [1] "s: "  "1534"
## [1] "s: "  "1535"
## [1] "s: "  "1536"
## [1] "s: "  "1537"
## [1] "s: "  "1538"
## [1] "s: "  "1539"
## [1] "s: "  "1540"
## [1] "s: "  "1541"
## [1] "s: "  "1542"
## [1] "s: "  "1543"
## [1] "s: "  "1544"
## [1] "s: "  "1545"
## [1] "s: "  "1546"
## [1] "s: "  "1547"
## [1] "s: "  "1548"
## [1] "s: "  "1549"
## [1] "s: "  "1550"
## [1] "s: "  "1551"
## [1] "s: "  "1552"
## [1] "s: "  "1553"
## [1] "s: "  "1554"
## [1] "s: "  "1555"
## [1] "s: "  "1556"
## [1] "s: "  "1557"
## [1] "s: "  "1558"
## [1] "s: "  "1559"
## [1] "s: "  "1560"
## [1] "s: "  "1561"
## [1] "s: "  "1562"
## [1] "s: "  "1563"
## [1] "s: "  "1564"
## [1] "s: "  "1565"
## [1] "s: "  "1566"
## [1] "s: "  "1567"
## [1] "s: "  "1568"
## [1] "s: "  "1569"
## [1] "s: "  "1570"
## [1] "s: "  "1571"
## [1] "s: "  "1572"
## [1] "s: "  "1573"
## [1] "s: "  "1574"
## [1] "s: "  "1575"
## [1] "s: "  "1576"
## [1] "s: "  "1577"
## [1] "s: "  "1578"
## [1] "s: "  "1579"
## [1] "s: "  "1580"
## [1] "s: "  "1581"
## [1] "s: "  "1582"
## [1] "s: "  "1583"
## [1] "s: "  "1584"
## [1] "s: "  "1585"
## [1] "s: "  "1586"
## [1] "s: "  "1587"
## [1] "s: "  "1588"
## [1] "s: "  "1589"
## [1] "s: "  "1590"
## [1] "s: "  "1591"
## [1] "s: "  "1592"
## [1] "s: "  "1593"
## [1] "s: "  "1594"
## [1] "s: "  "1595"
## [1] "s: "  "1596"
## [1] "s: "  "1597"
## [1] "s: "  "1598"
## [1] "s: "  "1599"
## [1] "s: "  "1600"
## [1] "s: "  "1601"
## [1] "s: "  "1602"
## [1] "s: "  "1603"
## [1] "s: "  "1604"
## [1] "s: "  "1605"
## [1] "s: "  "1606"
## [1] "s: "  "1607"
## [1] "s: "  "1608"
## [1] "s: "  "1609"
## [1] "s: "  "1610"
## [1] "s: "  "1611"
## [1] "s: "  "1612"
## [1] "s: "  "1613"
## [1] "s: "  "1614"
## [1] "s: "  "1615"
## [1] "s: "  "1616"
## [1] "s: "  "1617"
## [1] "s: "  "1618"
## [1] "s: "  "1619"
## [1] "s: "  "1620"
## [1] "s: "  "1621"
## [1] "s: "  "1622"
## [1] "s: "  "1623"
## [1] "s: "  "1624"
## [1] "s: "  "1625"
## [1] "s: "  "1626"
## [1] "s: "  "1627"
## [1] "s: "  "1628"
## [1] "s: "  "1629"
## [1] "s: "  "1630"
## [1] "s: "  "1631"
## [1] "s: "  "1632"
## [1] "s: "  "1633"
## [1] "s: "  "1634"
## [1] "s: "  "1635"
## [1] "s: "  "1636"
## [1] "s: "  "1637"
## [1] "s: "  "1638"
## [1] "s: "  "1639"
## [1] "s: "  "1640"
## [1] "s: "  "1641"
## [1] "s: "  "1642"
## [1] "s: "  "1643"
## [1] "s: "  "1644"
## [1] "s: "  "1645"
## [1] "s: "  "1646"
## [1] "s: "  "1647"
## [1] "s: "  "1648"
## [1] "s: "  "1649"
## [1] "s: "  "1650"
## [1] "s: "  "1651"
## [1] "s: "  "1652"
## [1] "s: "  "1653"
## [1] "s: "  "1654"
## [1] "s: "  "1655"
## [1] "s: "  "1656"
## [1] "s: "  "1657"
## [1] "s: "  "1658"
## [1] "s: "  "1659"
## [1] "s: "  "1660"
## [1] "s: "  "1661"
## [1] "s: "  "1662"
## [1] "s: "  "1663"
## [1] "s: "  "1664"
## [1] "s: "  "1665"
## [1] "s: "  "1666"
## [1] "s: "  "1667"
## [1] "s: "  "1668"
## [1] "s: "  "1669"
## [1] "s: "  "1670"
## [1] "s: "  "1671"
## [1] "s: "  "1672"
## [1] "s: "  "1673"
## [1] "s: "  "1674"
## [1] "s: "  "1675"
## [1] "s: "  "1676"
## [1] "s: "  "1677"
## [1] "s: "  "1678"
## [1] "s: "  "1679"
## [1] "s: "  "1680"
## [1] "s: "  "1681"
## [1] "s: "  "1682"
## [1] "s: "  "1683"
## [1] "s: "  "1684"
## [1] "s: "  "1685"
## [1] "s: "  "1686"
## [1] "s: "  "1687"
## [1] "s: "  "1688"
## [1] "s: "  "1689"
## [1] "s: "  "1690"
## [1] "s: "  "1691"
## [1] "s: "  "1692"
## [1] "s: "  "1693"
## [1] "s: "  "1694"
## [1] "s: "  "1695"
## [1] "s: "  "1696"
## [1] "s: "  "1697"
## [1] "s: "  "1698"
## [1] "s: "  "1699"
## [1] "s: "  "1700"
## [1] "s: "  "1701"
## [1] "s: "  "1702"
## [1] "s: "  "1703"
## [1] "s: "  "1704"
## [1] "s: "  "1705"
## [1] "s: "  "1706"
## [1] "s: "  "1707"
## [1] "s: "  "1708"
## [1] "s: "  "1709"
## [1] "s: "  "1710"
## [1] "s: "  "1711"
## [1] "s: "  "1712"
## [1] "s: "  "1713"
## [1] "s: "  "1714"
## [1] "s: "  "1715"
## [1] "s: "  "1716"
## [1] "s: "  "1717"
## [1] "s: "  "1718"
## [1] "s: "  "1719"
## [1] "s: "  "1720"
## [1] "s: "  "1721"
## [1] "s: "  "1722"
## [1] "s: "  "1723"
## [1] "s: "  "1724"
## [1] "s: "  "1725"
## [1] "s: "  "1726"
## [1] "s: "  "1727"
## [1] "s: "  "1728"
## [1] "s: "  "1729"
## [1] "s: "  "1730"
## [1] "s: "  "1731"
## [1] "s: "  "1732"
## [1] "s: "  "1733"
## [1] "s: "  "1734"
## [1] "s: "  "1735"
## [1] "s: "  "1736"
## [1] "s: "  "1737"
## [1] "s: "  "1738"
## [1] "s: "  "1739"
## [1] "s: "  "1740"
## [1] "s: "  "1741"
## [1] "s: "  "1742"
## [1] "s: "  "1743"
## [1] "s: "  "1744"
## [1] "s: "  "1745"
## [1] "s: "  "1746"
## [1] "s: "  "1747"
## [1] "s: "  "1748"
## [1] "s: "  "1749"
## [1] "s: "  "1750"
## [1] "s: "  "1751"
## [1] "s: "  "1752"
## [1] "s: "  "1753"
## [1] "s: "  "1754"
## [1] "s: "  "1755"
## [1] "s: "  "1756"
## [1] "s: "  "1757"
## [1] "s: "  "1758"
## [1] "s: "  "1759"
## [1] "s: "  "1760"
## [1] "s: "  "1761"
## [1] "s: "  "1762"
## [1] "s: "  "1763"
## [1] "s: "  "1764"
## [1] "s: "  "1765"
## [1] "s: "  "1766"
## [1] "s: "  "1767"
## [1] "s: "  "1768"
## [1] "s: "  "1769"
## [1] "s: "  "1770"
## [1] "s: "  "1771"
## [1] "s: "  "1772"
## [1] "s: "  "1773"
## [1] "s: "  "1774"
## [1] "s: "  "1775"
## [1] "s: "  "1776"
## [1] "s: "  "1777"
## [1] "s: "  "1778"
## [1] "s: "  "1779"
## [1] "s: "  "1780"
## [1] "s: "  "1781"
## [1] "s: "  "1782"
## [1] "s: "  "1783"
## [1] "s: "  "1784"
## [1] "s: "  "1785"
## [1] "s: "  "1786"
## [1] "s: "  "1787"
## [1] "s: "  "1788"
## [1] "s: "  "1789"
## [1] "s: "  "1790"
## [1] "s: "  "1791"
## [1] "s: "  "1792"
## [1] "s: "  "1793"
## [1] "s: "  "1794"
## [1] "s: "  "1795"
## [1] "s: "  "1796"
## [1] "s: "  "1797"
## [1] "s: "  "1798"
## [1] "s: "  "1799"
## [1] "s: "  "1800"
## [1] "s: "  "1801"
## [1] "s: "  "1802"
## [1] "s: "  "1803"
## [1] "s: "  "1804"
## [1] "s: "  "1805"
## [1] "s: "  "1806"
## [1] "s: "  "1807"
## [1] "s: "  "1808"
## [1] "s: "  "1809"
## [1] "s: "  "1810"
## [1] "s: "  "1811"
## [1] "s: "  "1812"
## [1] "s: "  "1813"
## [1] "s: "  "1814"
## [1] "s: "  "1815"
## [1] "s: "  "1816"
## [1] "s: "  "1817"
## [1] "s: "  "1818"
## [1] "s: "  "1819"
## [1] "s: "  "1820"
## [1] "s: "  "1821"
## [1] "s: "  "1822"
## [1] "s: "  "1823"
## [1] "s: "  "1824"
## [1] "s: "  "1825"
## [1] "s: "  "1826"
## [1] "s: "  "1827"
## [1] "s: "  "1828"
## [1] "s: "  "1829"
## [1] "s: "  "1830"
## [1] "s: "  "1831"
## [1] "s: "  "1832"
## [1] "s: "  "1833"
## [1] "s: "  "1834"
## [1] "s: "  "1835"
## [1] "s: "  "1836"
## [1] "s: "  "1837"
## [1] "s: "  "1838"
## [1] "s: "  "1839"
## [1] "s: "  "1840"
## [1] "s: "  "1841"
## [1] "s: "  "1842"
## [1] "s: "  "1843"
## [1] "s: "  "1844"
## [1] "s: "  "1845"
## [1] "s: "  "1846"
## [1] "s: "  "1847"
## [1] "s: "  "1848"
## [1] "s: "  "1849"
## [1] "s: "  "1850"
## [1] "s: "  "1851"
## [1] "s: "  "1852"
## [1] "s: "  "1853"
## [1] "s: "  "1854"
## [1] "s: "  "1855"
## [1] "s: "  "1856"
## [1] "s: "  "1857"
## [1] "s: "  "1858"
## [1] "s: "  "1859"
## [1] "s: "  "1860"
## [1] "s: "  "1861"
## [1] "s: "  "1862"
## [1] "s: "  "1863"
## [1] "s: "  "1864"
## [1] "s: "  "1865"
## [1] "s: "  "1866"
## [1] "s: "  "1867"
## [1] "s: "  "1868"
## [1] "s: "  "1869"
## [1] "s: "  "1870"
## [1] "s: "  "1871"
## [1] "s: "  "1872"
## [1] "s: "  "1873"
## [1] "s: "  "1874"
## [1] "s: "  "1875"
## [1] "s: "  "1876"
## [1] "s: "  "1877"
## [1] "s: "  "1878"
## [1] "s: "  "1879"
## [1] "s: "  "1880"
## [1] "s: "  "1881"
## [1] "s: "  "1882"
## [1] "s: "  "1883"
## [1] "s: "  "1884"
## [1] "s: "  "1885"
## [1] "s: "  "1886"
## [1] "s: "  "1887"
## [1] "s: "  "1888"
## [1] "s: "  "1889"
## [1] "s: "  "1890"
## [1] "s: "  "1891"
## [1] "s: "  "1892"
## [1] "s: "  "1893"
## [1] "s: "  "1894"
## [1] "s: "  "1895"
## [1] "s: "  "1896"
## [1] "s: "  "1897"
## [1] "s: "  "1898"
## [1] "s: "  "1899"
## [1] "s: "  "1900"
## [1] "s: "  "1901"
## [1] "s: "  "1902"
## [1] "s: "  "1903"
## [1] "s: "  "1904"
## [1] "s: "  "1905"
## [1] "s: "  "1906"
## [1] "s: "  "1907"
## [1] "s: "  "1908"
## [1] "s: "  "1909"
## [1] "s: "  "1910"
## [1] "s: "  "1911"
## [1] "s: "  "1912"
## [1] "s: "  "1913"
## [1] "s: "  "1914"
## [1] "s: "  "1915"
## [1] "s: "  "1916"
## [1] "s: "  "1917"
## [1] "s: "  "1918"
## [1] "s: "  "1919"
## [1] "s: "  "1920"
## [1] "s: "  "1921"
## [1] "s: "  "1922"
## [1] "s: "  "1923"
## [1] "s: "  "1924"
## [1] "s: "  "1925"
## [1] "s: "  "1926"
## [1] "s: "  "1927"
## [1] "s: "  "1928"
## [1] "s: "  "1929"
## [1] "s: "  "1930"
## [1] "s: "  "1931"
## [1] "s: "  "1932"
## [1] "s: "  "1933"
## [1] "s: "  "1934"
## [1] "s: "  "1935"
## [1] "s: "  "1936"
## [1] "s: "  "1937"
## [1] "s: "  "1938"
## [1] "s: "  "1939"
## [1] "s: "  "1940"
## [1] "s: "  "1941"
## [1] "s: "  "1942"
## [1] "s: "  "1943"
## [1] "s: "  "1944"
## [1] "s: "  "1945"
## [1] "s: "  "1946"
## [1] "s: "  "1947"
## [1] "s: "  "1948"
## [1] "s: "  "1949"
## [1] "s: "  "1950"
## [1] "s: "  "1951"
## [1] "s: "  "1952"
## [1] "s: "  "1953"
## [1] "s: "  "1954"
## [1] "s: "  "1955"
## [1] "s: "  "1956"
## [1] "s: "  "1957"
## [1] "s: "  "1958"
## [1] "s: "  "1959"
## [1] "s: "  "1960"
## [1] "s: "  "1961"
## [1] "s: "  "1962"
## [1] "s: "  "1963"
## [1] "s: "  "1964"
## [1] "s: "  "1965"
## [1] "s: "  "1966"
## [1] "s: "  "1967"
## [1] "s: "  "1968"
## [1] "s: "  "1969"
## [1] "s: "  "1970"
## [1] "s: "  "1971"
## [1] "s: "  "1972"
## [1] "s: "  "1973"
## [1] "s: "  "1974"
## [1] "s: "  "1975"
## [1] "s: "  "1976"
## [1] "s: "  "1977"
## [1] "s: "  "1978"
## [1] "s: "  "1979"
## [1] "s: "  "1980"
## [1] "s: "  "1981"
## [1] "s: "  "1982"
## [1] "s: "  "1983"
## [1] "s: "  "1984"
## [1] "s: "  "1985"
## [1] "s: "  "1986"
## [1] "s: "  "1987"
## [1] "s: "  "1988"
## [1] "s: "  "1989"
## [1] "s: "  "1990"
## [1] "s: "  "1991"
## [1] "s: "  "1992"
## [1] "s: "  "1993"
## [1] "s: "  "1994"
## [1] "s: "  "1995"
## [1] "s: "  "1996"
## [1] "s: "  "1997"
## [1] "s: "  "1998"
## [1] "s: "  "1999"
## [1] "s: "  "2000"
## [1] "s: "  "2001"
## [1] "s: "  "2002"
## [1] "s: "  "2003"
## [1] "s: "  "2004"
## [1] "s: "  "2005"
## [1] "s: "  "2006"
## [1] "s: "  "2007"
## [1] "s: "  "2008"
## [1] "s: "  "2009"
## [1] "s: "  "2010"
## [1] "s: "  "2011"
## [1] "s: "  "2012"
## [1] "s: "  "2013"
## [1] "s: "  "2014"
## [1] "s: "  "2015"
## [1] "s: "  "2016"
## [1] "s: "  "2017"
## [1] "s: "  "2018"
## [1] "s: "  "2019"
## [1] "s: "  "2020"
## [1] "s: "  "2021"
## [1] "s: "  "2022"
## [1] "s: "  "2023"
## [1] "s: "  "2024"
## [1] "s: "  "2025"
## [1] "s: "  "2026"
## [1] "s: "  "2027"
## [1] "s: "  "2028"
## [1] "s: "  "2029"
## [1] "s: "  "2030"
## [1] "s: "  "2031"
## [1] "s: "  "2032"
## [1] "s: "  "2033"
## [1] "s: "  "2034"
## [1] "s: "  "2035"
## [1] "s: "  "2036"
## [1] "s: "  "2037"
## [1] "s: "  "2038"
## [1] "s: "  "2039"
## [1] "s: "  "2040"
## [1] "s: "  "2041"
## [1] "s: "  "2042"
## [1] "s: "  "2043"
## [1] "s: "  "2044"
## [1] "s: "  "2045"
## [1] "s: "  "2046"
## [1] "s: "  "2047"
## [1] "s: "  "2048"
## [1] "s: "  "2049"
## [1] "s: "  "2050"
## [1] "s: "  "2051"
## [1] "s: "  "2052"
## [1] "s: "  "2053"
## [1] "s: "  "2054"
## [1] "s: "  "2055"
## [1] "s: "  "2056"
## [1] "s: "  "2057"
## [1] "s: "  "2058"
## [1] "s: "  "2059"
## [1] "s: "  "2060"
## [1] "s: "  "2061"
## [1] "s: "  "2062"
## [1] "s: "  "2063"
## [1] "s: "  "2064"
## [1] "s: "  "2065"
## [1] "s: "  "2066"
## [1] "s: "  "2067"
## [1] "s: "  "2068"
## [1] "s: "  "2069"
## [1] "s: "  "2070"
## [1] "s: "  "2071"
## [1] "s: "  "2072"
## [1] "s: "  "2073"
## [1] "s: "  "2074"
## [1] "s: "  "2075"
## [1] "s: "  "2076"
## [1] "s: "  "2077"
## [1] "s: "  "2078"
## [1] "s: "  "2079"
## [1] "s: "  "2080"
## [1] "s: "  "2081"
## [1] "s: "  "2082"
## [1] "s: "  "2083"
## [1] "s: "  "2084"
## [1] "s: "  "2085"
## [1] "s: "  "2086"
## [1] "s: "  "2087"
## [1] "s: "  "2088"
## [1] "s: "  "2089"
## [1] "s: "  "2090"
## [1] "s: "  "2091"
## [1] "s: "  "2092"
## [1] "s: "  "2093"
## [1] "s: "  "2094"
## [1] "s: "  "2095"
## [1] "s: "  "2096"
## [1] "s: "  "2097"
## [1] "s: "  "2098"
## [1] "s: "  "2099"
## [1] "s: "  "2100"
## [1] "s: "  "2101"
## [1] "s: "  "2102"
## [1] "s: "  "2103"
## [1] "s: "  "2104"
## [1] "s: "  "2105"
## [1] "s: "  "2106"
## [1] "s: "  "2107"
## [1] "s: "  "2108"
## [1] "s: "  "2109"
## [1] "s: "  "2110"
## [1] "s: "  "2111"
## [1] "s: "  "2112"
## [1] "s: "  "2113"
## [1] "s: "  "2114"
## [1] "s: "  "2115"
## [1] "s: "  "2116"
## [1] "s: "  "2117"
## [1] "s: "  "2118"
## [1] "s: "  "2119"
## [1] "s: "  "2120"
## [1] "s: "  "2121"
## [1] "s: "  "2122"
## [1] "s: "  "2123"
## [1] "s: "  "2124"
## [1] "s: "  "2125"
## [1] "s: "  "2126"
## [1] "s: "  "2127"
## [1] "s: "  "2128"
## [1] "s: "  "2129"
## [1] "s: "  "2130"
## [1] "s: "  "2131"
## [1] "s: "  "2132"
## [1] "s: "  "2133"
## [1] "s: "  "2134"
## [1] "s: "  "2135"
## [1] "s: "  "2136"
## [1] "s: "  "2137"
## [1] "s: "  "2138"
## [1] "s: "  "2139"
## [1] "s: "  "2140"
## [1] "s: "  "2141"
## [1] "s: "  "2142"
## [1] "s: "  "2143"
## [1] "s: "  "2144"
## [1] "s: "  "2145"
## [1] "s: "  "2146"
## [1] "s: "  "2147"
## [1] "s: "  "2148"
## [1] "s: "  "2149"
## [1] "s: "  "2150"
## [1] "s: "  "2151"
## [1] "s: "  "2152"
## [1] "s: "  "2153"
## [1] "s: "  "2154"
## [1] "s: "  "2155"
## [1] "s: "  "2156"
## [1] "s: "  "2157"
## [1] "s: "  "2158"
## [1] "s: "  "2159"
## [1] "s: "  "2160"
## [1] "s: "  "2161"
## [1] "s: "  "2162"
## [1] "s: "  "2163"
## [1] "s: "  "2164"
## [1] "s: "  "2165"
## [1] "s: "  "2166"
## [1] "s: "  "2167"
## [1] "s: "  "2168"
## [1] "s: "  "2169"
## [1] "s: "  "2170"
## [1] "s: "  "2171"
## [1] "s: "  "2172"
## [1] "s: "  "2173"
## [1] "s: "  "2174"
## [1] "s: "  "2175"
## [1] "s: "  "2176"
## [1] "s: "  "2177"
## [1] "s: "  "2178"
## [1] "s: "  "2179"
## [1] "s: "  "2180"
## [1] "s: "  "2181"
## [1] "s: "  "2182"
## [1] "s: "  "2183"
## [1] "s: "  "2184"
## [1] "s: "  "2185"
## [1] "s: "  "2186"
## [1] "s: "  "2187"
## [1] "s: "  "2188"
## [1] "s: "  "2189"
## [1] "s: "  "2190"
## [1] "s: "  "2191"
## [1] "s: "  "2192"
## [1] "s: "  "2193"
## [1] "s: "  "2194"
## [1] "s: "  "2195"
## [1] "s: "  "2196"
## [1] "s: "  "2197"
## [1] "s: "  "2198"
## [1] "s: "  "2199"
## [1] "s: "  "2200"
## [1] "s: "  "2201"
## [1] "s: "  "2202"
## [1] "s: "  "2203"
## [1] "s: "  "2204"
## [1] "s: "  "2205"
## [1] "s: "  "2206"
## [1] "s: "  "2207"
## [1] "s: "  "2208"
## [1] "s: "  "2209"
## [1] "s: "  "2210"
## [1] "s: "  "2211"
## [1] "s: "  "2212"
## [1] "s: "  "2213"
## [1] "s: "  "2214"
## [1] "s: "  "2215"
## [1] "s: "  "2216"
## [1] "s: "  "2217"
## [1] "s: "  "2218"
## [1] "s: "  "2219"
## [1] "s: "  "2220"
## [1] "s: "  "2221"
## [1] "s: "  "2222"
## [1] "s: "  "2223"
## [1] "s: "  "2224"
## [1] "s: "  "2225"
## [1] "s: "  "2226"
## [1] "s: "  "2227"
## [1] "s: "  "2228"
## [1] "s: "  "2229"
## [1] "s: "  "2230"
## [1] "s: "  "2231"
## [1] "s: "  "2232"
## [1] "s: "  "2233"
## [1] "s: "  "2234"
## [1] "s: "  "2235"
## [1] "s: "  "2236"
## [1] "s: "  "2237"
## [1] "s: "  "2238"
## [1] "s: "  "2239"
## [1] "s: "  "2240"
## [1] "s: "  "2241"
## [1] "s: "  "2242"
## [1] "s: "  "2243"
## [1] "s: "  "2244"
## [1] "s: "  "2245"
## [1] "s: "  "2246"
## [1] "s: "  "2247"
## [1] "s: "  "2248"
## [1] "s: "  "2249"
## [1] "s: "  "2250"
## [1] "s: "  "2251"
## [1] "s: "  "2252"
## [1] "s: "  "2253"
## [1] "s: "  "2254"
## [1] "s: "  "2255"
## [1] "s: "  "2256"
## [1] "s: "  "2257"
## [1] "s: "  "2258"
## [1] "s: "  "2259"
## [1] "s: "  "2260"
## [1] "s: "  "2261"
## [1] "s: "  "2262"
## [1] "s: "  "2263"
## [1] "s: "  "2264"
## [1] "s: "  "2265"
## [1] "s: "  "2266"
## [1] "s: "  "2267"
## [1] "s: "  "2268"
## [1] "s: "  "2269"
## [1] "s: "  "2270"
## [1] "s: "  "2271"
## [1] "s: "  "2272"
## [1] "s: "  "2273"
## [1] "s: "  "2274"
## [1] "s: "  "2275"
## [1] "s: "  "2276"
## [1] "s: "  "2277"
## [1] "s: "  "2278"
## [1] "s: "  "2279"
## [1] "s: "  "2280"
## [1] "s: "  "2281"
## [1] "s: "  "2282"
## [1] "s: "  "2283"
## [1] "s: "  "2284"
## [1] "s: "  "2285"
## [1] "s: "  "2286"
## [1] "s: "  "2287"
## [1] "s: "  "2288"
## [1] "s: "  "2289"
## [1] "s: "  "2290"
## [1] "s: "  "2291"
## [1] "s: "  "2292"
## [1] "s: "  "2293"
## [1] "s: "  "2294"
## [1] "s: "  "2295"
## [1] "s: "  "2296"
## [1] "s: "  "2297"
## [1] "s: "  "2298"
## [1] "s: "  "2299"
## [1] "s: "  "2300"
## [1] "s: "  "2301"
## [1] "s: "  "2302"
## [1] "s: "  "2303"
## [1] "s: "  "2304"
## [1] "s: "  "2305"
## [1] "s: "  "2306"
## [1] "s: "  "2307"
## [1] "s: "  "2308"
## [1] "s: "  "2309"
## [1] "s: "  "2310"
## [1] "s: "  "2311"
## [1] "s: "  "2312"
## [1] "s: "  "2313"
## [1] "s: "  "2314"
## [1] "s: "  "2315"
## [1] "s: "  "2316"
## [1] "s: "  "2317"
## [1] "s: "  "2318"
## [1] "s: "  "2319"
## [1] "s: "  "2320"
## [1] "s: "  "2321"
## [1] "s: "  "2322"
## [1] "s: "  "2323"
## [1] "s: "  "2324"
## [1] "s: "  "2325"
## [1] "s: "  "2326"
## [1] "s: "  "2327"
## [1] "s: "  "2328"
## [1] "s: "  "2329"
## [1] "s: "  "2330"
## [1] "s: "  "2331"
## [1] "s: "  "2332"
## [1] "s: "  "2333"
## [1] "s: "  "2334"
## [1] "s: "  "2335"
## [1] "s: "  "2336"
## [1] "s: "  "2337"
## [1] "s: "  "2338"
## [1] "s: "  "2339"
## [1] "s: "  "2340"
## [1] "s: "  "2341"
## [1] "s: "  "2342"
## [1] "s: "  "2343"
## [1] "s: "  "2344"
## [1] "s: "  "2345"
## [1] "s: "  "2346"
## [1] "s: "  "2347"
## [1] "s: "  "2348"
## [1] "s: "  "2349"
## [1] "s: "  "2350"
## [1] "s: "  "2351"
## [1] "s: "  "2352"
## [1] "s: "  "2353"
## [1] "s: "  "2354"
## [1] "s: "  "2355"
## [1] "s: "  "2356"
## [1] "s: "  "2357"
## [1] "s: "  "2358"
## [1] "s: "  "2359"
## [1] "s: "  "2360"
## [1] "s: "  "2361"
## [1] "s: "  "2362"
## [1] "s: "  "2363"
## [1] "s: "  "2364"
## [1] "s: "  "2365"
## [1] "s: "  "2366"
## [1] "s: "  "2367"
## [1] "s: "  "2368"
## [1] "s: "  "2369"
## [1] "s: "  "2370"
## [1] "s: "  "2371"
## [1] "s: "  "2372"
## [1] "s: "  "2373"
## [1] "s: "  "2374"
## [1] "s: "  "2375"
## [1] "s: "  "2376"
## [1] "s: "  "2377"
## [1] "s: "  "2378"
## [1] "s: "  "2379"
## [1] "s: "  "2380"
## [1] "s: "  "2381"
## [1] "s: "  "2382"
## [1] "s: "  "2383"
## [1] "s: "  "2384"
## [1] "s: "  "2385"
## [1] "s: "  "2386"
## [1] "s: "  "2387"
## [1] "s: "  "2388"
## [1] "s: "  "2389"
## [1] "s: "  "2390"
## [1] "s: "  "2391"
## [1] "s: "  "2392"
## [1] "s: "  "2393"
## [1] "s: "  "2394"
## [1] "s: "  "2395"
## [1] "s: "  "2396"
## [1] "s: "  "2397"
## [1] "s: "  "2398"
## [1] "s: "  "2399"
## [1] "s: "  "2400"
## [1] "s: "  "2401"
## [1] "s: "  "2402"
## [1] "s: "  "2403"
## [1] "s: "  "2404"
## [1] "s: "  "2405"
## [1] "s: "  "2406"
## [1] "s: "  "2407"
## [1] "s: "  "2408"
## [1] "s: "  "2409"
## [1] "s: "  "2410"
## [1] "s: "  "2411"
## [1] "s: "  "2412"
## [1] "s: "  "2413"
## [1] "s: "  "2414"
## [1] "s: "  "2415"
## [1] "s: "  "2416"
## [1] "s: "  "2417"
## [1] "s: "  "2418"
## [1] "s: "  "2419"
## [1] "s: "  "2420"
## [1] "s: "  "2421"
## [1] "s: "  "2422"
## [1] "s: "  "2423"
## [1] "s: "  "2424"
## [1] "s: "  "2425"
## [1] "s: "  "2426"
## [1] "s: "  "2427"
## [1] "s: "  "2428"
## [1] "s: "  "2429"
## [1] "s: "  "2430"
## [1] "s: "  "2431"
## [1] "s: "  "2432"
## [1] "s: "  "2433"
## [1] "s: "  "2434"
## [1] "s: "  "2435"
## [1] "s: "  "2436"
## [1] "s: "  "2437"
## [1] "s: "  "2438"
## [1] "s: "  "2439"
## [1] "s: "  "2440"
## [1] "s: "  "2441"
## [1] "s: "  "2442"
## [1] "s: "  "2443"
## [1] "s: "  "2444"
## [1] "s: "  "2445"
## [1] "s: "  "2446"
## [1] "s: "  "2447"
## [1] "s: "  "2448"
## [1] "s: "  "2449"
## [1] "s: "  "2450"
## [1] "s: "  "2451"
## [1] "s: "  "2452"
## [1] "s: "  "2453"
## [1] "s: "  "2454"
## [1] "s: "  "2455"
## [1] "s: "  "2456"
## [1] "s: "  "2457"
## [1] "s: "  "2458"
## [1] "s: "  "2459"
## [1] "s: "  "2460"
## [1] "s: "  "2461"
## [1] "s: "  "2462"
## [1] "s: "  "2463"
## [1] "s: "  "2464"
## [1] "s: "  "2465"
## [1] "s: "  "2466"
## [1] "s: "  "2467"
## [1] "s: "  "2468"
## [1] "s: "  "2469"
## [1] "s: "  "2470"
## [1] "s: "  "2471"
## [1] "s: "  "2472"
## [1] "s: "  "2473"
## [1] "s: "  "2474"
## [1] "s: "  "2475"
## [1] "s: "  "2476"
## [1] "s: "  "2477"
## [1] "s: "  "2478"
## [1] "s: "  "2479"
## [1] "s: "  "2480"
## [1] "s: "  "2481"
## [1] "s: "  "2482"
## [1] "s: "  "2483"
## [1] "s: "  "2484"
## [1] "s: "  "2485"
## [1] "s: "  "2486"
## [1] "s: "  "2487"
## [1] "s: "  "2488"
## [1] "s: "  "2489"
## [1] "s: "  "2490"
## [1] "s: "  "2491"
## [1] "s: "  "2492"
## [1] "s: "  "2493"
## [1] "s: "  "2494"
## [1] "s: "  "2495"
## [1] "s: "  "2496"
## [1] "s: "  "2497"
## [1] "s: "  "2498"
## [1] "s: "  "2499"
## [1] "s: "  "2500"
## [1] "s: "  "2501"
## [1] "s: "  "2502"
## [1] "s: "  "2503"
## [1] "s: "  "2504"
## [1] "s: "  "2505"
## [1] "s: "  "2506"
## [1] "s: "  "2507"
## [1] "s: "  "2508"
## [1] "s: "  "2509"
## [1] "s: "  "2510"
## [1] "s: "  "2511"
## [1] "s: "  "2512"
## [1] "s: "  "2513"
## [1] "s: "  "2514"
## [1] "s: "  "2515"
## [1] "s: "  "2516"
## [1] "s: "  "2517"
## [1] "s: "  "2518"
## [1] "s: "  "2519"
## [1] "s: "  "2520"
## [1] "s: "  "2521"
## [1] "s: "  "2522"
## [1] "s: "  "2523"
## [1] "s: "  "2524"
## [1] "s: "  "2525"
## [1] "s: "  "2526"
## [1] "s: "  "2527"
## [1] "s: "  "2528"
## [1] "s: "  "2529"
## [1] "s: "  "2530"
## [1] "s: "  "2531"
## [1] "s: "  "2532"
## [1] "s: "  "2533"
## [1] "s: "  "2534"
## [1] "s: "  "2535"
## [1] "s: "  "2536"
## [1] "s: "  "2537"
## [1] "s: "  "2538"
## [1] "s: "  "2539"
## [1] "s: "  "2540"
## [1] "s: "  "2541"
## [1] "s: "  "2542"
## [1] "s: "  "2543"
## [1] "s: "  "2544"
## [1] "s: "  "2545"
## [1] "s: "  "2546"
## [1] "s: "  "2547"
## [1] "s: "  "2548"
## [1] "s: "  "2549"
## [1] "s: "  "2550"
## [1] "s: "  "2551"
## [1] "s: "  "2552"
## [1] "s: "  "2553"
## [1] "s: "  "2554"
## [1] "s: "  "2555"
## [1] "s: "  "2556"
## [1] "s: "  "2557"
## [1] "s: "  "2558"
## [1] "s: "  "2559"
## [1] "s: "  "2560"
## [1] "s: "  "2561"
## [1] "s: "  "2562"
## [1] "s: "  "2563"
## [1] "s: "  "2564"
## [1] "s: "  "2565"
## [1] "s: "  "2566"
## [1] "s: "  "2567"
## [1] "s: "  "2568"
## [1] "s: "  "2569"
## [1] "s: "  "2570"
## [1] "s: "  "2571"
## [1] "s: "  "2572"
## [1] "s: "  "2573"
## [1] "s: "  "2574"
## [1] "s: "  "2575"
## [1] "s: "  "2576"
## [1] "s: "  "2577"
## [1] "s: "  "2578"
## [1] "s: "  "2579"
## [1] "s: "  "2580"
## [1] "s: "  "2581"
## [1] "s: "  "2582"
## [1] "s: "  "2583"
## [1] "s: "  "2584"
## [1] "s: "  "2585"
## [1] "s: "  "2586"
## [1] "s: "  "2587"
## [1] "s: "  "2588"
## [1] "s: "  "2589"
## [1] "s: "  "2590"
## [1] "s: "  "2591"
## [1] "s: "  "2592"
## [1] "s: "  "2593"
## [1] "s: "  "2594"
## [1] "s: "  "2595"
## [1] "s: "  "2596"
## [1] "s: "  "2597"
## [1] "s: "  "2598"
## [1] "s: "  "2599"
## [1] "s: "  "2600"
## [1] "s: "  "2601"
## [1] "s: "  "2602"
## [1] "s: "  "2603"
## [1] "s: "  "2604"
## [1] "s: "  "2605"
## [1] "s: "  "2606"
## [1] "s: "  "2607"
## [1] "s: "  "2608"
## [1] "s: "  "2609"
## [1] "s: "  "2610"
## [1] "s: "  "2611"
## [1] "s: "  "2612"
## [1] "s: "  "2613"
## [1] "s: "  "2614"
## [1] "s: "  "2615"
## [1] "s: "  "2616"
## [1] "s: "  "2617"
## [1] "s: "  "2618"
## [1] "s: "  "2619"
## [1] "s: "  "2620"
## [1] "s: "  "2621"
## [1] "s: "  "2622"
## [1] "s: "  "2623"
## [1] "s: "  "2624"
## [1] "s: "  "2625"
## [1] "s: "  "2626"
## [1] "s: "  "2627"
## [1] "s: "  "2628"
## [1] "s: "  "2629"
## [1] "s: "  "2630"
## [1] "s: "  "2631"
## [1] "s: "  "2632"
## [1] "s: "  "2633"
## [1] "s: "  "2634"
## [1] "s: "  "2635"
## [1] "s: "  "2636"
## [1] "s: "  "2637"
## [1] "s: "  "2638"
## [1] "s: "  "2639"
## [1] "s: "  "2640"
## [1] "s: "  "2641"
## [1] "s: "  "2642"
## [1] "s: "  "2643"
## [1] "s: "  "2644"
## [1] "s: "  "2645"
## [1] "s: "  "2646"
## [1] "s: "  "2647"
## [1] "s: "  "2648"
## [1] "s: "  "2649"
## [1] "s: "  "2650"
## [1] "s: "  "2651"
## [1] "s: "  "2652"
## [1] "s: "  "2653"
## [1] "s: "  "2654"
## [1] "s: "  "2655"
## [1] "s: "  "2656"
## [1] "s: "  "2657"
## [1] "s: "  "2658"
## [1] "s: "  "2659"
## [1] "s: "  "2660"
## [1] "s: "  "2661"
## [1] "s: "  "2662"
## [1] "s: "  "2663"
## [1] "s: "  "2664"
## [1] "s: "  "2665"
## [1] "s: "  "2666"
## [1] "s: "  "2667"
## [1] "s: "  "2668"
## [1] "s: "  "2669"
## [1] "s: "  "2670"
## [1] "s: "  "2671"
## [1] "s: "  "2672"
## [1] "s: "  "2673"
## [1] "s: "  "2674"
## [1] "s: "  "2675"
## [1] "s: "  "2676"
## [1] "s: "  "2677"
## [1] "s: "  "2678"
## [1] "s: "  "2679"
## [1] "s: "  "2680"
## [1] "s: "  "2681"
## [1] "s: "  "2682"
## [1] "s: "  "2683"
## [1] "s: "  "2684"
## [1] "s: "  "2685"
## [1] "s: "  "2686"
## [1] "s: "  "2687"
## [1] "s: "  "2688"
## [1] "s: "  "2689"
## [1] "s: "  "2690"
## [1] "s: "  "2691"
## [1] "s: "  "2692"
## [1] "s: "  "2693"
## [1] "s: "  "2694"
## [1] "s: "  "2695"
## [1] "s: "  "2696"
## [1] "s: "  "2697"
## [1] "s: "  "2698"
## [1] "s: "  "2699"
## [1] "s: "  "2700"
## [1] "s: "  "2701"
## [1] "s: "  "2702"
## [1] "s: "  "2703"
## [1] "s: "  "2704"
## [1] "s: "  "2705"
## [1] "s: "  "2706"
## [1] "s: "  "2707"
## [1] "s: "  "2708"
## [1] "s: "  "2709"
## [1] "s: "  "2710"
## [1] "s: "  "2711"
## [1] "s: "  "2712"
## [1] "s: "  "2713"
## [1] "s: "  "2714"
## [1] "s: "  "2715"
## [1] "s: "  "2716"
## [1] "s: "  "2717"
## [1] "s: "  "2718"
## [1] "s: "  "2719"
## [1] "s: "  "2720"
## [1] "s: "  "2721"
## [1] "s: "  "2722"
## [1] "s: "  "2723"
## [1] "s: "  "2724"
## [1] "s: "  "2725"
## [1] "s: "  "2726"
## [1] "s: "  "2727"
## [1] "s: "  "2728"
## [1] "s: "  "2729"
## [1] "s: "  "2730"
## [1] "s: "  "2731"
## [1] "s: "  "2732"
## [1] "s: "  "2733"
## [1] "s: "  "2734"
## [1] "s: "  "2735"
## [1] "s: "  "2736"
## [1] "s: "  "2737"
## [1] "s: "  "2738"
## [1] "s: "  "2739"
## [1] "s: "  "2740"
## [1] "s: "  "2741"
## [1] "s: "  "2742"
## [1] "s: "  "2743"
## [1] "s: "  "2744"
## [1] "s: "  "2745"
## [1] "s: "  "2746"
## [1] "s: "  "2747"
## [1] "s: "  "2748"
## [1] "s: "  "2749"
## [1] "s: "  "2750"
## [1] "s: "  "2751"
## [1] "s: "  "2752"
## [1] "s: "  "2753"
## [1] "s: "  "2754"
## [1] "s: "  "2755"
## [1] "s: "  "2756"
## [1] "s: "  "2757"
## [1] "s: "  "2758"
## [1] "s: "  "2759"
## [1] "s: "  "2760"
## [1] "s: "  "2761"
## [1] "s: "  "2762"
## [1] "s: "  "2763"
## [1] "s: "  "2764"
## [1] "s: "  "2765"
## [1] "s: "  "2766"
## [1] "s: "  "2767"
## [1] "s: "  "2768"
## [1] "s: "  "2769"
## [1] "s: "  "2770"
## [1] "s: "  "2771"
## [1] "s: "  "2772"
## [1] "s: "  "2773"
## [1] "s: "  "2774"
## [1] "s: "  "2775"
## [1] "s: "  "2776"
## [1] "s: "  "2777"
## [1] "s: "  "2778"
## [1] "s: "  "2779"
## [1] "s: "  "2780"
## [1] "s: "  "2781"
## [1] "s: "  "2782"
## [1] "s: "  "2783"
## [1] "s: "  "2784"
## [1] "s: "  "2785"
## [1] "s: "  "2786"
## [1] "s: "  "2787"
## [1] "s: "  "2788"
## [1] "s: "  "2789"
## [1] "s: "  "2790"
## [1] "s: "  "2791"
## [1] "s: "  "2792"
## [1] "s: "  "2793"
## [1] "s: "  "2794"
## [1] "s: "  "2795"
## [1] "s: "  "2796"
## [1] "s: "  "2797"
## [1] "s: "  "2798"
## [1] "s: "  "2799"
## [1] "s: "  "2800"
## [1] "s: "  "2801"
## [1] "s: "  "2802"
## [1] "s: "  "2803"
## [1] "s: "  "2804"
## [1] "s: "  "2805"
## [1] "s: "  "2806"
## [1] "s: "  "2807"
## [1] "s: "  "2808"
## [1] "s: "  "2809"
## [1] "s: "  "2810"
## [1] "s: "  "2811"
## [1] "s: "  "2812"
## [1] "s: "  "2813"
## [1] "s: "  "2814"
## [1] "s: "  "2815"
## [1] "s: "  "2816"
## [1] "s: "  "2817"
## [1] "s: "  "2818"
## [1] "s: "  "2819"
## [1] "s: "  "2820"
## [1] "s: "  "2821"
## [1] "s: "  "2822"
## [1] "s: "  "2823"
## [1] "s: "  "2824"
## [1] "s: "  "2825"
## [1] "s: "  "2826"
## [1] "s: "  "2827"
## [1] "s: "  "2828"
## [1] "s: "  "2829"
## [1] "s: "  "2830"
## [1] "s: "  "2831"
## [1] "s: "  "2832"
## [1] "s: "  "2833"
## [1] "s: "  "2834"
## [1] "s: "  "2835"
## [1] "s: "  "2836"
## [1] "s: "  "2837"
## [1] "s: "  "2838"
## [1] "s: "  "2839"
## [1] "s: "  "2840"
## [1] "s: "  "2841"
## [1] "s: "  "2842"
## [1] "s: "  "2843"
## [1] "s: "  "2844"
## [1] "s: "  "2845"
## [1] "s: "  "2846"
## [1] "s: "  "2847"
## [1] "s: "  "2848"
## [1] "s: "  "2849"
## [1] "s: "  "2850"
## [1] "s: "  "2851"
## [1] "s: "  "2852"
## [1] "s: "  "2853"
## [1] "s: "  "2854"
## [1] "s: "  "2855"
## [1] "s: "  "2856"
## [1] "s: "  "2857"
## [1] "s: "  "2858"
## [1] "s: "  "2859"
## [1] "s: "  "2860"
## [1] "s: "  "2861"
## [1] "s: "  "2862"
## [1] "s: "  "2863"
## [1] "s: "  "2864"
## [1] "s: "  "2865"
## [1] "s: "  "2866"
## [1] "s: "  "2867"
## [1] "s: "  "2868"
## [1] "s: "  "2869"
## [1] "s: "  "2870"
## [1] "s: "  "2871"
## [1] "s: "  "2872"
## [1] "s: "  "2873"
## [1] "s: "  "2874"
## [1] "s: "  "2875"
## [1] "s: "  "2876"
## [1] "s: "  "2877"
## [1] "s: "  "2878"
## [1] "s: "  "2879"
## [1] "s: "  "2880"
## [1] "s: "  "2881"
## [1] "s: "  "2882"
## [1] "s: "  "2883"
## [1] "s: "  "2884"
## [1] "s: "  "2885"
## [1] "s: "  "2886"
## [1] "s: "  "2887"
## [1] "s: "  "2888"
## [1] "s: "  "2889"
## [1] "s: "  "2890"
## [1] "s: "  "2891"
## [1] "s: "  "2892"
## [1] "s: "  "2893"
## [1] "s: "  "2894"
## [1] "s: "  "2895"
## [1] "s: "  "2896"
## [1] "s: "  "2897"
## [1] "s: "  "2898"
## [1] "s: "  "2899"
## [1] "s: "  "2900"
## [1] "s: "  "2901"
## [1] "s: "  "2902"
## [1] "s: "  "2903"
## [1] "s: "  "2904"
## [1] "s: "  "2905"
## [1] "s: "  "2906"
## [1] "s: "  "2907"
## [1] "s: "  "2908"
## [1] "s: "  "2909"
## [1] "s: "  "2910"
## [1] "s: "  "2911"
## [1] "s: "  "2912"
## [1] "s: "  "2913"
## [1] "s: "  "2914"
## [1] "s: "  "2915"
## [1] "s: "  "2916"
## [1] "s: "  "2917"
## [1] "s: "  "2918"
## [1] "s: "  "2919"
## [1] "s: "  "2920"
## [1] "s: "  "2921"
## [1] "s: "  "2922"
## [1] "s: "  "2923"
## [1] "s: "  "2924"
## [1] "s: "  "2925"
## [1] "s: "  "2926"
## [1] "s: "  "2927"
## [1] "s: "  "2928"
## [1] "s: "  "2929"
## [1] "s: "  "2930"
## [1] "s: "  "2931"
## [1] "s: "  "2932"
## [1] "s: "  "2933"
## [1] "s: "  "2934"
## [1] "s: "  "2935"
## [1] "s: "  "2936"
## [1] "s: "  "2937"
## [1] "s: "  "2938"
## [1] "s: "  "2939"
## [1] "s: "  "2940"
## [1] "s: "  "2941"
## [1] "s: "  "2942"
## [1] "s: "  "2943"
## [1] "s: "  "2944"
## [1] "s: "  "2945"
## [1] "s: "  "2946"
## [1] "s: "  "2947"
## [1] "s: "  "2948"
## [1] "s: "  "2949"
## [1] "s: "  "2950"
## [1] "s: "  "2951"
## [1] "s: "  "2952"
## [1] "s: "  "2953"
## [1] "s: "  "2954"
## [1] "s: "  "2955"
## [1] "s: "  "2956"
## [1] "s: "  "2957"
## [1] "s: "  "2958"
## [1] "s: "  "2959"
## [1] "s: "  "2960"
## [1] "s: "  "2961"
## [1] "s: "  "2962"
## [1] "s: "  "2963"
## [1] "s: "  "2964"
## [1] "s: "  "2965"
## [1] "s: "  "2966"
## [1] "s: "  "2967"
## [1] "s: "  "2968"
## [1] "s: "  "2969"
## [1] "s: "  "2970"
## [1] "s: "  "2971"
## [1] "s: "  "2972"
## [1] "s: "  "2973"
## [1] "s: "  "2974"
## [1] "s: "  "2975"
## [1] "s: "  "2976"
## [1] "s: "  "2977"
## [1] "s: "  "2978"
## [1] "s: "  "2979"
## [1] "s: "  "2980"
## [1] "s: "  "2981"
## [1] "s: "  "2982"
## [1] "s: "  "2983"
## [1] "s: "  "2984"
## [1] "s: "  "2985"
## [1] "s: "  "2986"
## [1] "s: "  "2987"
## [1] "s: "  "2988"
## [1] "s: "  "2989"
## [1] "s: "  "2990"
## [1] "s: "  "2991"
## [1] "s: "  "2992"
## [1] "s: "  "2993"
## [1] "s: "  "2994"
## [1] "s: "  "2995"
## [1] "s: "  "2996"
## [1] "s: "  "2997"
## [1] "s: "  "2998"
## [1] "s: "  "2999"
## [1] "s: "  "3000"
## [1] "s: "  "3001"
## [1] "s: "  "3002"
## [1] "s: "  "3003"
## [1] "s: "  "3004"
## [1] "s: "  "3005"
## [1] "s: "  "3006"
## [1] "s: "  "3007"
## [1] "s: "  "3008"
## [1] "s: "  "3009"
## [1] "s: "  "3010"
## [1] "s: "  "3011"
## [1] "s: "  "3012"
## [1] "s: "  "3013"
## [1] "s: "  "3014"
## [1] "s: "  "3015"
## [1] "s: "  "3016"
## [1] "s: "  "3017"
## [1] "s: "  "3018"
## [1] "s: "  "3019"
## [1] "s: "  "3020"
## [1] "s: "  "3021"
## [1] "s: "  "3022"
## [1] "s: "  "3023"
## [1] "s: "  "3024"
## [1] "s: "  "3025"
## [1] "s: "  "3026"
## [1] "s: "  "3027"
## [1] "s: "  "3028"
## [1] "s: "  "3029"
## [1] "s: "  "3030"
## [1] "s: "  "3031"
## [1] "s: "  "3032"
## [1] "s: "  "3033"
## [1] "s: "  "3034"
## [1] "s: "  "3035"
## [1] "s: "  "3036"
## [1] "s: "  "3037"
## [1] "s: "  "3038"
## [1] "s: "  "3039"
## [1] "s: "  "3040"
## [1] "s: "  "3041"
## [1] "s: "  "3042"
## [1] "s: "  "3043"
## [1] "s: "  "3044"
## [1] "s: "  "3045"
## [1] "s: "  "3046"
## [1] "s: "  "3047"
## [1] "s: "  "3048"
## [1] "s: "  "3049"
## [1] "s: "  "3050"
## [1] "s: "  "3051"
## [1] "s: "  "3052"
## [1] "s: "  "3053"
## [1] "s: "  "3054"
## [1] "s: "  "3055"
## [1] "s: "  "3056"
## [1] "s: "  "3057"
## [1] "s: "  "3058"
## [1] "s: "  "3059"
## [1] "s: "  "3060"
## [1] "s: "  "3061"
## [1] "s: "  "3062"
## [1] "s: "  "3063"
## [1] "s: "  "3064"
## [1] "s: "  "3065"
## [1] "s: "  "3066"
## [1] "s: "  "3067"
## [1] "s: "  "3068"
## [1] "s: "  "3069"
## [1] "s: "  "3070"
## [1] "s: "  "3071"
## [1] "s: "  "3072"
## [1] "s: "  "3073"
## [1] "s: "  "3074"
## [1] "s: "  "3075"
## [1] "s: "  "3076"
## [1] "s: "  "3077"
## [1] "s: "  "3078"
## [1] "s: "  "3079"
## [1] "s: "  "3080"
## [1] "s: "  "3081"
## [1] "s: "  "3082"
## [1] "s: "  "3083"
## [1] "s: "  "3084"
## [1] "s: "  "3085"
## [1] "s: "  "3086"
## [1] "s: "  "3087"
## [1] "s: "  "3088"
## [1] "s: "  "3089"
## [1] "s: "  "3090"
## [1] "s: "  "3091"
## [1] "s: "  "3092"
## [1] "s: "  "3093"
## [1] "s: "  "3094"
## [1] "s: "  "3095"
## [1] "s: "  "3096"
## [1] "s: "  "3097"
## [1] "s: "  "3098"
## [1] "s: "  "3099"
## [1] "s: "  "3100"
## [1] "s: "  "3101"
## [1] "s: "  "3102"
## [1] "s: "  "3103"
## [1] "s: "  "3104"
## [1] "s: "  "3105"
## [1] "s: "  "3106"
## [1] "s: "  "3107"
## [1] "s: "  "3108"
## [1] "s: "  "3109"
## [1] "s: "  "3110"
## [1] "s: "  "3111"
## [1] "s: "  "3112"
## [1] "s: "  "3113"
## [1] "s: "  "3114"
## [1] "s: "  "3115"
## [1] "s: "  "3116"
## [1] "s: "  "3117"
## [1] "s: "  "3118"
## [1] "s: "  "3119"
## [1] "s: "  "3120"
## [1] "s: "  "3121"
## [1] "s: "  "3122"
## [1] "s: "  "3123"
## [1] "s: "  "3124"
## [1] "s: "  "3125"
## [1] "s: "  "3126"
## [1] "s: "  "3127"
## [1] "s: "  "3128"
## [1] "s: "  "3129"
## [1] "s: "  "3130"
## [1] "s: "  "3131"
## [1] "s: "  "3132"
## [1] "s: "  "3133"
## [1] "s: "  "3134"
## [1] "s: "  "3135"
## [1] "s: "  "3136"
## [1] "s: "  "3137"
## [1] "s: "  "3138"
## [1] "s: "  "3139"
## [1] "s: "  "3140"
## [1] "s: "  "3141"
## [1] "s: "  "3142"
## [1] "s: "  "3143"
## [1] "s: "  "3144"
## [1] "s: "  "3145"
## [1] "s: "  "3146"
## [1] "s: "  "3147"
## [1] "s: "  "3148"
## [1] "s: "  "3149"
## [1] "s: "  "3150"
## [1] "s: "  "3151"
## [1] "s: "  "3152"
## [1] "s: "  "3153"
## [1] "s: "  "3154"
## [1] "s: "  "3155"
## [1] "s: "  "3156"
## [1] "s: "  "3157"
## [1] "s: "  "3158"
## [1] "s: "  "3159"
## [1] "s: "  "3160"
## [1] "s: "  "3161"
## [1] "s: "  "3162"
## [1] "s: "  "3163"
## [1] "s: "  "3164"
## [1] "s: "  "3165"
## [1] "s: "  "3166"
## [1] "s: "  "3167"
## [1] "s: "  "3168"
## [1] "s: "  "3169"
## [1] "s: "  "3170"
## [1] "s: "  "3171"
## [1] "s: "  "3172"
## [1] "s: "  "3173"
## [1] "s: "  "3174"
## [1] "s: "  "3175"
## [1] "s: "  "3176"
## [1] "s: "  "3177"
## [1] "s: "  "3178"
## [1] "s: "  "3179"
## [1] "s: "  "3180"
## [1] "s: "  "3181"
## [1] "s: "  "3182"
## [1] "s: "  "3183"
## [1] "s: "  "3184"
## [1] "s: "  "3185"
## [1] "s: "  "3186"
## [1] "s: "  "3187"
## [1] "s: "  "3188"
## [1] "s: "  "3189"
## [1] "s: "  "3190"
## [1] "s: "  "3191"
## [1] "s: "  "3192"
## [1] "s: "  "3193"
## [1] "s: "  "3194"
## [1] "s: "  "3195"
## [1] "s: "  "3196"
## [1] "s: "  "3197"
## [1] "s: "  "3198"
## [1] "s: "  "3199"
## [1] "s: "  "3200"
## [1] "s: "  "3201"
## [1] "s: "  "3202"
## [1] "s: "  "3203"
## [1] "s: "  "3204"
## [1] "s: "  "3205"
## [1] "s: "  "3206"
## [1] "s: "  "3207"
## [1] "s: "  "3208"
## [1] "s: "  "3209"
## [1] "s: "  "3210"
## [1] "s: "  "3211"
## [1] "s: "  "3212"
## [1] "s: "  "3213"
## [1] "s: "  "3214"
## [1] "s: "  "3215"
## [1] "s: "  "3216"
## [1] "s: "  "3217"
## [1] "s: "  "3218"
## [1] "s: "  "3219"
## [1] "s: "  "3220"
## [1] "s: "  "3221"
## [1] "s: "  "3222"
## [1] "s: "  "3223"
## [1] "s: "  "3224"
## [1] "s: "  "3225"
## [1] "s: "  "3226"
## [1] "s: "  "3227"
## [1] "s: "  "3228"
## [1] "s: "  "3229"
## [1] "s: "  "3230"
## [1] "s: "  "3231"
## [1] "s: "  "3232"
## [1] "s: "  "3233"
## [1] "s: "  "3234"
## [1] "s: "  "3235"
## [1] "s: "  "3236"
## [1] "s: "  "3237"
## [1] "s: "  "3238"
## [1] "s: "  "3239"
## [1] "s: "  "3240"
## [1] "s: "  "3241"
## [1] "s: "  "3242"
## [1] "s: "  "3243"
## [1] "s: "  "3244"
## [1] "s: "  "3245"
## [1] "s: "  "3246"
## [1] "s: "  "3247"
## [1] "s: "  "3248"
## [1] "s: "  "3249"
## [1] "s: "  "3250"
## [1] "s: "  "3251"
## [1] "s: "  "3252"
## [1] "s: "  "3253"
## [1] "s: "  "3254"
## [1] "s: "  "3255"
## [1] "s: "  "3256"
## [1] "s: "  "3257"
## [1] "s: "  "3258"
## [1] "s: "  "3259"
## [1] "s: "  "3260"
## [1] "s: "  "3261"
## [1] "s: "  "3262"
## [1] "s: "  "3263"
## [1] "s: "  "3264"
## [1] "s: "  "3265"
## [1] "s: "  "3266"
## [1] "s: "  "3267"
## [1] "s: "  "3268"
## [1] "s: "  "3269"
## [1] "s: "  "3270"
## [1] "s: "  "3271"
## [1] "s: "  "3272"
## [1] "s: "  "3273"
## [1] "s: "  "3274"
## [1] "s: "  "3275"
## [1] "s: "  "3276"
## [1] "s: "  "3277"
## [1] "s: "  "3278"
## [1] "s: "  "3279"
## [1] "s: "  "3280"
## [1] "s: "  "3281"
## [1] "s: "  "3282"
## [1] "s: "  "3283"
## [1] "s: "  "3284"
## [1] "s: "  "3285"
## [1] "s: "  "3286"
## [1] "s: "  "3287"
## [1] "s: "  "3288"
## [1] "s: "  "3289"
## [1] "s: "  "3290"
## [1] "s: "  "3291"
## [1] "s: "  "3292"
## [1] "s: "  "3293"
## [1] "s: "  "3294"
## [1] "s: "  "3295"
## [1] "s: "  "3296"
## [1] "s: "  "3297"
## [1] "s: "  "3298"
## [1] "s: "  "3299"
## [1] "s: "  "3300"
## [1] "s: "  "3301"
## [1] "s: "  "3302"
## [1] "s: "  "3303"
## [1] "s: "  "3304"
## [1] "s: "  "3305"
## [1] "s: "  "3306"
## [1] "s: "  "3307"
## [1] "s: "  "3308"
## [1] "s: "  "3309"
## [1] "s: "  "3310"
## [1] "s: "  "3311"
## [1] "s: "  "3312"
## [1] "s: "  "3313"
## [1] "s: "  "3314"
## [1] "s: "  "3315"
## [1] "s: "  "3316"
## [1] "s: "  "3317"
## [1] "s: "  "3318"
## [1] "s: "  "3319"
## [1] "s: "  "3320"
## [1] "s: "  "3321"
## [1] "s: "  "3322"
## [1] "s: "  "3323"
## [1] "s: "  "3324"
## [1] "s: "  "3325"
## [1] "s: "  "3326"
## [1] "s: "  "3327"
## [1] "s: "  "3328"
## [1] "s: "  "3329"
## [1] "s: "  "3330"
## [1] "s: "  "3331"
## [1] "s: "  "3332"
## [1] "s: "  "3333"
## [1] "s: "  "3334"
## [1] "s: "  "3335"
## [1] "s: "  "3336"
## [1] "s: "  "3337"
## [1] "s: "  "3338"
## [1] "s: "  "3339"
## [1] "s: "  "3340"
## [1] "s: "  "3341"
## [1] "s: "  "3342"
## [1] "s: "  "3343"
## [1] "s: "  "3344"
## [1] "s: "  "3345"
## [1] "s: "  "3346"
## [1] "s: "  "3347"
## [1] "s: "  "3348"
## [1] "s: "  "3349"
## [1] "s: "  "3350"
## [1] "s: "  "3351"
## [1] "s: "  "3352"
## [1] "s: "  "3353"
## [1] "s: "  "3354"
## [1] "s: "  "3355"
## [1] "s: "  "3356"
## [1] "s: "  "3357"
## [1] "s: "  "3358"
## [1] "s: "  "3359"
## [1] "s: "  "3360"
## [1] "s: "  "3361"
## [1] "s: "  "3362"
## [1] "s: "  "3363"
## [1] "s: "  "3364"
## [1] "s: "  "3365"
## [1] "s: "  "3366"
## [1] "s: "  "3367"
## [1] "s: "  "3368"
## [1] "s: "  "3369"
## [1] "s: "  "3370"
## [1] "s: "  "3371"
## [1] "s: "  "3372"
## [1] "s: "  "3373"
## [1] "s: "  "3374"
## [1] "s: "  "3375"
## [1] "s: "  "3376"
## [1] "s: "  "3377"
## [1] "s: "  "3378"
## [1] "s: "  "3379"
## [1] "s: "  "3380"
## [1] "s: "  "3381"
## [1] "s: "  "3382"
## [1] "s: "  "3383"
## [1] "s: "  "3384"
## [1] "s: "  "3385"
## [1] "s: "  "3386"
## [1] "s: "  "3387"
## [1] "s: "  "3388"
## [1] "s: "  "3389"
## [1] "s: "  "3390"
## [1] "s: "  "3391"
## [1] "s: "  "3392"
## [1] "s: "  "3393"
## [1] "s: "  "3394"
## [1] "s: "  "3395"
## [1] "s: "  "3396"
## [1] "s: "  "3397"
## [1] "s: "  "3398"
## [1] "s: "  "3399"
## [1] "s: "  "3400"
## [1] "s: "  "3401"
## [1] "s: "  "3402"
## [1] "s: "  "3403"
## [1] "s: "  "3404"
## [1] "s: "  "3405"
## [1] "s: "  "3406"
## [1] "s: "  "3407"
## [1] "s: "  "3408"
## [1] "s: "  "3409"
## [1] "s: "  "3410"
## [1] "s: "  "3411"
## [1] "s: "  "3412"
## [1] "s: "  "3413"
## [1] "s: "  "3414"
## [1] "s: "  "3415"
## [1] "s: "  "3416"
## [1] "s: "  "3417"
## [1] "s: "  "3418"
## [1] "s: "  "3419"
## [1] "s: "  "3420"
## [1] "s: "  "3421"
## [1] "s: "  "3422"
## [1] "s: "  "3423"
## [1] "s: "  "3424"
## [1] "s: "  "3425"
## [1] "s: "  "3426"
## [1] "s: "  "3427"
## [1] "s: "  "3428"
## [1] "s: "  "3429"
## [1] "s: "  "3430"
## [1] "s: "  "3431"
## [1] "s: "  "3432"
## [1] "s: "  "3433"
## [1] "s: "  "3434"
## [1] "s: "  "3435"
## [1] "s: "  "3436"
## [1] "s: "  "3437"
## [1] "s: "  "3438"
## [1] "s: "  "3439"
## [1] "s: "  "3440"
## [1] "s: "  "3441"
## [1] "s: "  "3442"
## [1] "s: "  "3443"
## [1] "s: "  "3444"
## [1] "s: "  "3445"
## [1] "s: "  "3446"
## [1] "s: "  "3447"
## [1] "s: "  "3448"
## [1] "s: "  "3449"
## [1] "s: "  "3450"
## [1] "s: "  "3451"
## [1] "s: "  "3452"
## [1] "s: "  "3453"
## [1] "s: "  "3454"
## [1] "s: "  "3455"
## [1] "s: "  "3456"
## [1] "s: "  "3457"
## [1] "s: "  "3458"
## [1] "s: "  "3459"
## [1] "s: "  "3460"
## [1] "s: "  "3461"
## [1] "s: "  "3462"
## [1] "s: "  "3463"
## [1] "s: "  "3464"
## [1] "s: "  "3465"
## [1] "s: "  "3466"
## [1] "s: "  "3467"
## [1] "s: "  "3468"
## [1] "s: "  "3469"
## [1] "s: "  "3470"
## [1] "s: "  "3471"
## [1] "s: "  "3472"
## [1] "s: "  "3473"
## [1] "s: "  "3474"
## [1] "s: "  "3475"
## [1] "s: "  "3476"
## [1] "s: "  "3477"
## [1] "s: "  "3478"
## [1] "s: "  "3479"
## [1] "s: "  "3480"
## [1] "s: "  "3481"
## [1] "s: "  "3482"
## [1] "s: "  "3483"
## [1] "s: "  "3484"
## [1] "s: "  "3485"
## [1] "s: "  "3486"
## [1] "s: "  "3487"
## [1] "s: "  "3488"
## [1] "s: "  "3489"
## [1] "s: "  "3490"
## [1] "s: "  "3491"
## [1] "s: "  "3492"
## [1] "s: "  "3493"
## [1] "s: "  "3494"
## [1] "s: "  "3495"
## [1] "s: "  "3496"
## [1] "s: "  "3497"
## [1] "s: "  "3498"
## [1] "s: "  "3499"
## [1] "s: "  "3500"
## [1] "s: "  "3501"
## [1] "s: "  "3502"
## [1] "s: "  "3503"
## [1] "s: "  "3504"
## [1] "s: "  "3505"
## [1] "s: "  "3506"
## [1] "s: "  "3507"
## [1] "s: "  "3508"
## [1] "s: "  "3509"
## [1] "s: "  "3510"
## [1] "s: "  "3511"
## [1] "s: "  "3512"
## [1] "s: "  "3513"
## [1] "s: "  "3514"
## [1] "s: "  "3515"
## [1] "s: "  "3516"
## [1] "s: "  "3517"
## [1] "s: "  "3518"
## [1] "s: "  "3519"
## [1] "s: "  "3520"
## [1] "s: "  "3521"
## [1] "s: "  "3522"
## [1] "s: "  "3523"
## [1] "s: "  "3524"
## [1] "s: "  "3525"
## [1] "s: "  "3526"
## [1] "s: "  "3527"
## [1] "s: "  "3528"
## [1] "s: "  "3529"
## [1] "s: "  "3530"
## [1] "s: "  "3531"
## [1] "s: "  "3532"
## [1] "s: "  "3533"
## [1] "s: "  "3534"
## [1] "s: "  "3535"
## [1] "s: "  "3536"
## [1] "s: "  "3537"
## [1] "s: "  "3538"
## [1] "s: "  "3539"
## [1] "s: "  "3540"
## [1] "s: "  "3541"
## [1] "s: "  "3542"
## [1] "s: "  "3543"
## [1] "s: "  "3544"
## [1] "s: "  "3545"
## [1] "s: "  "3546"
## [1] "s: "  "3547"
## [1] "s: "  "3548"
## [1] "s: "  "3549"
## [1] "s: "  "3550"
## [1] "s: "  "3551"
## [1] "s: "  "3552"
## [1] "s: "  "3553"
## [1] "s: "  "3554"
## [1] "s: "  "3555"
## [1] "s: "  "3556"
## [1] "s: "  "3557"
## [1] "s: "  "3558"
## [1] "s: "  "3559"
## [1] "s: "  "3560"
## [1] "s: "  "3561"
## [1] "s: "  "3562"
## [1] "s: "  "3563"
## [1] "s: "  "3564"
## [1] "s: "  "3565"
## [1] "s: "  "3566"
## [1] "s: "  "3567"
## [1] "s: "  "3568"
## [1] "s: "  "3569"
## [1] "s: "  "3570"
## [1] "s: "  "3571"
## [1] "s: "  "3572"
## [1] "s: "  "3573"
## [1] "s: "  "3574"
## [1] "s: "  "3575"
## [1] "s: "  "3576"
## [1] "s: "  "3577"
## [1] "s: "  "3578"
## [1] "s: "  "3579"
## [1] "s: "  "3580"
## [1] "s: "  "3581"
## [1] "s: "  "3582"
## [1] "s: "  "3583"
## [1] "s: "  "3584"
## [1] "s: "  "3585"
## [1] "s: "  "3586"
## [1] "s: "  "3587"
## [1] "s: "  "3588"
## [1] "s: "  "3589"
## [1] "s: "  "3590"
## [1] "s: "  "3591"
## [1] "s: "  "3592"
## [1] "s: "  "3593"
## [1] "s: "  "3594"
## [1] "s: "  "3595"
## [1] "s: "  "3596"
## [1] "s: "  "3597"
## [1] "s: "  "3598"
## [1] "s: "  "3599"
## [1] "s: "  "3600"
## [1] "s: "  "3601"
## [1] "s: "  "3602"
## [1] "s: "  "3603"
## [1] "s: "  "3604"
## [1] "s: "  "3605"
## [1] "s: "  "3606"
## [1] "s: "  "3607"
## [1] "s: "  "3608"
## [1] "s: "  "3609"
## [1] "s: "  "3610"
## [1] "s: "  "3611"
## [1] "s: "  "3612"
## [1] "s: "  "3613"
## [1] "s: "  "3614"
## [1] "s: "  "3615"
## [1] "s: "  "3616"
## [1] "s: "  "3617"
## [1] "s: "  "3618"
## [1] "s: "  "3619"
## [1] "s: "  "3620"
## [1] "s: "  "3621"
## [1] "s: "  "3622"
## [1] "s: "  "3623"
## [1] "s: "  "3624"
## [1] "s: "  "3625"
## [1] "s: "  "3626"
## [1] "s: "  "3627"
## [1] "s: "  "3628"
## [1] "s: "  "3629"
## [1] "s: "  "3630"
## [1] "s: "  "3631"
## [1] "s: "  "3632"
## [1] "s: "  "3633"
## [1] "s: "  "3634"
## [1] "s: "  "3635"
## [1] "s: "  "3636"
## [1] "s: "  "3637"
## [1] "s: "  "3638"
## [1] "s: "  "3639"
## [1] "s: "  "3640"
## [1] "s: "  "3641"
## [1] "s: "  "3642"
## [1] "s: "  "3643"
## [1] "s: "  "3644"
## [1] "s: "  "3645"
## [1] "s: "  "3646"
## [1] "s: "  "3647"
## [1] "s: "  "3648"
## [1] "s: "  "3649"
## [1] "s: "  "3650"
## [1] "s: "  "3651"
## [1] "s: "  "3652"
## [1] "s: "  "3653"
## [1] "s: "  "3654"
## [1] "s: "  "3655"
## [1] "s: "  "3656"
## [1] "s: "  "3657"
## [1] "s: "  "3658"
## [1] "s: "  "3659"
## [1] "s: "  "3660"
## [1] "s: "  "3661"
## [1] "s: "  "3662"
## [1] "s: "  "3663"
## [1] "s: "  "3664"
## [1] "s: "  "3665"
## [1] "s: "  "3666"
## [1] "s: "  "3667"
## [1] "s: "  "3668"
## [1] "s: "  "3669"
## [1] "s: "  "3670"
## [1] "s: "  "3671"
## [1] "s: "  "3672"
## [1] "s: "  "3673"
## [1] "s: "  "3674"
## [1] "s: "  "3675"
## [1] "s: "  "3676"
## [1] "s: "  "3677"
## [1] "s: "  "3678"
## [1] "s: "  "3679"
## [1] "s: "  "3680"
## [1] "s: "  "3681"
## [1] "s: "  "3682"
## [1] "s: "  "3683"
## [1] "s: "  "3684"
## [1] "s: "  "3685"
## [1] "s: "  "3686"
## [1] "s: "  "3687"
## [1] "s: "  "3688"
## [1] "s: "  "3689"
## [1] "s: "  "3690"
## [1] "s: "  "3691"
## [1] "s: "  "3692"
## [1] "s: "  "3693"
## [1] "s: "  "3694"
## [1] "s: "  "3695"
## [1] "s: "  "3696"
## [1] "s: "  "3697"
## [1] "s: "  "3698"
## [1] "s: "  "3699"
## [1] "s: "  "3700"
## [1] "s: "  "3701"
## [1] "s: "  "3702"
## [1] "s: "  "3703"
## [1] "s: "  "3704"
## [1] "s: "  "3705"
## [1] "s: "  "3706"
## [1] "s: "  "3707"
## [1] "s: "  "3708"
## [1] "s: "  "3709"
## [1] "s: "  "3710"
## [1] "s: "  "3711"
## [1] "s: "  "3712"
## [1] "s: "  "3713"
## [1] "s: "  "3714"
## [1] "s: "  "3715"
## [1] "s: "  "3716"
## [1] "s: "  "3717"
## [1] "s: "  "3718"
## [1] "s: "  "3719"
## [1] "s: "  "3720"
## [1] "s: "  "3721"
## [1] "s: "  "3722"
## [1] "s: "  "3723"
## [1] "s: "  "3724"
## [1] "s: "  "3725"
## [1] "s: "  "3726"
## [1] "s: "  "3727"
## [1] "s: "  "3728"
## [1] "s: "  "3729"
## [1] "s: "  "3730"
## [1] "s: "  "3731"
## [1] "s: "  "3732"
## [1] "s: "  "3733"
## [1] "s: "  "3734"
## [1] "s: "  "3735"
## [1] "s: "  "3736"
## [1] "s: "  "3737"
## [1] "s: "  "3738"
## [1] "s: "  "3739"
## [1] "s: "  "3740"
## [1] "s: "  "3741"
## [1] "s: "  "3742"
## [1] "s: "  "3743"
## [1] "s: "  "3744"
## [1] "s: "  "3745"
## [1] "s: "  "3746"
## [1] "s: "  "3747"
## [1] "s: "  "3748"
## [1] "s: "  "3749"
## [1] "s: "  "3750"
## [1] "s: "  "3751"
## [1] "s: "  "3752"
## [1] "s: "  "3753"
## [1] "s: "  "3754"
## [1] "s: "  "3755"
## [1] "s: "  "3756"
## [1] "s: "  "3757"
## [1] "s: "  "3758"
## [1] "s: "  "3759"
## [1] "s: "  "3760"
## [1] "s: "  "3761"
## [1] "s: "  "3762"
## [1] "s: "  "3763"
## [1] "s: "  "3764"
## [1] "s: "  "3765"
## [1] "s: "  "3766"
## [1] "s: "  "3767"
## [1] "s: "  "3768"
## [1] "s: "  "3769"
## [1] "s: "  "3770"
## [1] "s: "  "3771"
## [1] "s: "  "3772"
## [1] "s: "  "3773"
## [1] "s: "  "3774"
## [1] "s: "  "3775"
## [1] "s: "  "3776"
## [1] "s: "  "3777"
## [1] "s: "  "3778"
## [1] "s: "  "3779"
## [1] "s: "  "3780"
## [1] "s: "  "3781"
## [1] "s: "  "3782"
## [1] "s: "  "3783"
## [1] "s: "  "3784"
## [1] "s: "  "3785"
## [1] "s: "  "3786"
## [1] "s: "  "3787"
## [1] "s: "  "3788"
## [1] "s: "  "3789"
## [1] "s: "  "3790"
## [1] "s: "  "3791"
## [1] "s: "  "3792"
## [1] "s: "  "3793"
## [1] "s: "  "3794"
## [1] "s: "  "3795"
## [1] "s: "  "3796"
## [1] "s: "  "3797"
## [1] "s: "  "3798"
## [1] "s: "  "3799"
## [1] "s: "  "3800"
## [1] "s: "  "3801"
## [1] "s: "  "3802"
## [1] "s: "  "3803"
## [1] "s: "  "3804"
## [1] "s: "  "3805"
## [1] "s: "  "3806"
## [1] "s: "  "3807"
## [1] "s: "  "3808"
## [1] "s: "  "3809"
## [1] "s: "  "3810"
## [1] "s: "  "3811"
## [1] "s: "  "3812"
## [1] "s: "  "3813"
## [1] "s: "  "3814"
## [1] "s: "  "3815"
## [1] "s: "  "3816"
## [1] "s: "  "3817"
## [1] "s: "  "3818"
## [1] "s: "  "3819"
## [1] "s: "  "3820"
## [1] "s: "  "3821"
## [1] "s: "  "3822"
## [1] "s: "  "3823"
## [1] "s: "  "3824"
## [1] "s: "  "3825"
## [1] "s: "  "3826"
## [1] "s: "  "3827"
## [1] "s: "  "3828"
## [1] "s: "  "3829"
## [1] "s: "  "3830"
## [1] "s: "  "3831"
## [1] "s: "  "3832"
## [1] "s: "  "3833"
## [1] "s: "  "3834"
## [1] "s: "  "3835"
## [1] "s: "  "3836"
## [1] "s: "  "3837"
## [1] "s: "  "3838"
## [1] "s: "  "3839"
## [1] "s: "  "3840"
## [1] "s: "  "3841"
## [1] "s: "  "3842"
## [1] "s: "  "3843"
## [1] "s: "  "3844"
## [1] "s: "  "3845"
## [1] "s: "  "3846"
## [1] "s: "  "3847"
## [1] "s: "  "3848"
## [1] "s: "  "3849"
## [1] "s: "  "3850"
## [1] "s: "  "3851"
## [1] "s: "  "3852"
## [1] "s: "  "3853"
## [1] "s: "  "3854"
## [1] "s: "  "3855"
## [1] "s: "  "3856"
## [1] "s: "  "3857"
## [1] "s: "  "3858"
## [1] "s: "  "3859"
## [1] "s: "  "3860"
## [1] "s: "  "3861"
## [1] "s: "  "3862"
## [1] "s: "  "3863"
## [1] "s: "  "3864"
## [1] "s: "  "3865"
## [1] "s: "  "3866"
## [1] "s: "  "3867"
## [1] "s: "  "3868"
## [1] "s: "  "3869"
## [1] "s: "  "3870"
## [1] "s: "  "3871"
## [1] "s: "  "3872"
## [1] "s: "  "3873"
## [1] "s: "  "3874"
## [1] "s: "  "3875"
## [1] "s: "  "3876"
## [1] "s: "  "3877"
## [1] "s: "  "3878"
## [1] "s: "  "3879"
## [1] "s: "  "3880"
## [1] "s: "  "3881"
## [1] "s: "  "3882"
## [1] "s: "  "3883"
## [1] "s: "  "3884"
## [1] "s: "  "3885"
## [1] "s: "  "3886"
## [1] "s: "  "3887"
## [1] "s: "  "3888"
## [1] "s: "  "3889"
## [1] "s: "  "3890"
## [1] "s: "  "3891"
## [1] "s: "  "3892"
## [1] "s: "  "3893"
## [1] "s: "  "3894"
## [1] "s: "  "3895"
## [1] "s: "  "3896"
## [1] "s: "  "3897"
## [1] "s: "  "3898"
## [1] "s: "  "3899"
## [1] "s: "  "3900"
## [1] "s: "  "3901"
## [1] "s: "  "3902"
## [1] "s: "  "3903"
## [1] "s: "  "3904"
## [1] "s: "  "3905"
## [1] "s: "  "3906"
## [1] "s: "  "3907"
## [1] "s: "  "3908"
## [1] "s: "  "3909"
## [1] "s: "  "3910"
## [1] "s: "  "3911"
## [1] "s: "  "3912"
## [1] "s: "  "3913"
## [1] "s: "  "3914"
## [1] "s: "  "3915"
## [1] "s: "  "3916"
## [1] "s: "  "3917"
## [1] "s: "  "3918"
## [1] "s: "  "3919"
## [1] "s: "  "3920"
## [1] "s: "  "3921"
## [1] "s: "  "3922"
## [1] "s: "  "3923"
## [1] "s: "  "3924"
## [1] "s: "  "3925"
## [1] "s: "  "3926"
## [1] "s: "  "3927"
## [1] "s: "  "3928"
## [1] "s: "  "3929"
## [1] "s: "  "3930"
## [1] "s: "  "3931"
## [1] "s: "  "3932"
## [1] "s: "  "3933"
## [1] "s: "  "3934"
## [1] "s: "  "3935"
## [1] "s: "  "3936"
## [1] "s: "  "3937"
## [1] "s: "  "3938"
## [1] "s: "  "3939"
## [1] "s: "  "3940"
## [1] "s: "  "3941"
## [1] "s: "  "3942"
## [1] "s: "  "3943"
## [1] "s: "  "3944"
## [1] "s: "  "3945"
## [1] "s: "  "3946"
## [1] "s: "  "3947"
## [1] "s: "  "3948"
## [1] "s: "  "3949"
## [1] "s: "  "3950"
## [1] "s: "  "3951"
## [1] "s: "  "3952"
## [1] "s: "  "3953"
## [1] "s: "  "3954"
## [1] "s: "  "3955"
## [1] "s: "  "3956"
## [1] "s: "  "3957"
## [1] "s: "  "3958"
## [1] "s: "  "3959"
## [1] "s: "  "3960"
## [1] "s: "  "3961"
## [1] "s: "  "3962"
## [1] "s: "  "3963"
## [1] "s: "  "3964"
## [1] "s: "  "3965"
## [1] "s: "  "3966"
## [1] "s: "  "3967"
## [1] "s: "  "3968"
## [1] "s: "  "3969"
## [1] "s: "  "3970"
## [1] "s: "  "3971"
## [1] "s: "  "3972"
## [1] "s: "  "3973"
## [1] "s: "  "3974"
## [1] "s: "  "3975"
## [1] "s: "  "3976"
## [1] "s: "  "3977"
## [1] "s: "  "3978"
## [1] "s: "  "3979"
## [1] "s: "  "3980"
## [1] "s: "  "3981"
## [1] "s: "  "3982"
## [1] "s: "  "3983"
## [1] "s: "  "3984"
## [1] "s: "  "3985"
## [1] "s: "  "3986"
## [1] "s: "  "3987"
## [1] "s: "  "3988"
## [1] "s: "  "3989"
## [1] "s: "  "3990"
## [1] "s: "  "3991"
## [1] "s: "  "3992"
## [1] "s: "  "3993"
## [1] "s: "  "3994"
## [1] "s: "  "3995"
## [1] "s: "  "3996"
## [1] "s: "  "3997"
## [1] "s: "  "3998"
## [1] "s: "  "3999"
## [1] "s: "  "4000"
## [1] "s: "  "4001"
## [1] "s: "  "4002"
## [1] "s: "  "4003"
## [1] "s: "  "4004"
## [1] "s: "  "4005"
## [1] "s: "  "4006"
## [1] "s: "  "4007"
## [1] "s: "  "4008"
## [1] "s: "  "4009"
## [1] "s: "  "4010"
## [1] "s: "  "4011"
## [1] "s: "  "4012"
## [1] "s: "  "4013"
## [1] "s: "  "4014"
## [1] "s: "  "4015"
## [1] "s: "  "4016"
## [1] "s: "  "4017"
## [1] "s: "  "4018"
## [1] "s: "  "4019"
## [1] "s: "  "4020"
## [1] "s: "  "4021"
## [1] "s: "  "4022"
## [1] "s: "  "4023"
## [1] "s: "  "4024"
## [1] "s: "  "4025"
## [1] "s: "  "4026"
## [1] "s: "  "4027"
## [1] "s: "  "4028"
## [1] "s: "  "4029"
## [1] "s: "  "4030"
## [1] "s: "  "4031"
## [1] "s: "  "4032"
## [1] "s: "  "4033"
## [1] "s: "  "4034"
## [1] "s: "  "4035"
## [1] "s: "  "4036"
## [1] "s: "  "4037"
## [1] "s: "  "4038"
## [1] "s: "  "4039"
## [1] "s: "  "4040"
## [1] "s: "  "4041"
## [1] "s: "  "4042"
## [1] "s: "  "4043"
## [1] "s: "  "4044"
## [1] "s: "  "4045"
## [1] "s: "  "4046"
## [1] "s: "  "4047"
## [1] "s: "  "4048"
## [1] "s: "  "4049"
## [1] "s: "  "4050"
## [1] "s: "  "4051"
## [1] "s: "  "4052"
## [1] "s: "  "4053"
## [1] "s: "  "4054"
## [1] "s: "  "4055"
## [1] "s: "  "4056"
## [1] "s: "  "4057"
## [1] "s: "  "4058"
## [1] "s: "  "4059"
## [1] "s: "  "4060"
## [1] "s: "  "4061"
## [1] "s: "  "4062"
## [1] "s: "  "4063"
## [1] "s: "  "4064"
## [1] "s: "  "4065"
## [1] "s: "  "4066"
## [1] "s: "  "4067"
## [1] "s: "  "4068"
## [1] "s: "  "4069"
## [1] "s: "  "4070"
## [1] "s: "  "4071"
## [1] "s: "  "4072"
## [1] "s: "  "4073"
## [1] "s: "  "4074"
## [1] "s: "  "4075"
## [1] "s: "  "4076"
## [1] "s: "  "4077"
## [1] "s: "  "4078"
## [1] "s: "  "4079"
## [1] "s: "  "4080"
## [1] "s: "  "4081"
## [1] "s: "  "4082"
## [1] "s: "  "4083"
## [1] "s: "  "4084"
## [1] "s: "  "4085"
## [1] "s: "  "4086"
## [1] "s: "  "4087"
## [1] "s: "  "4088"
## [1] "s: "  "4089"
## [1] "s: "  "4090"
## [1] "s: "  "4091"
## [1] "s: "  "4092"
## [1] "s: "  "4093"
## [1] "s: "  "4094"
## [1] "s: "  "4095"
## [1] "s: "  "4096"
## [1] "s: "  "4097"
## [1] "s: "  "4098"
## [1] "s: "  "4099"
## [1] "s: "  "4100"
## [1] "s: "  "4101"
## [1] "s: "  "4102"
## [1] "s: "  "4103"
## [1] "s: "  "4104"
## [1] "s: "  "4105"
## [1] "s: "  "4106"
## [1] "s: "  "4107"
## [1] "s: "  "4108"
## [1] "s: "  "4109"
## [1] "s: "  "4110"
## [1] "s: "  "4111"
## [1] "s: "  "4112"
## [1] "s: "  "4113"
## [1] "s: "  "4114"
## [1] "s: "  "4115"
## [1] "s: "  "4116"
## [1] "s: "  "4117"
## [1] "s: "  "4118"
## [1] "s: "  "4119"
## [1] "s: "  "4120"
## [1] "s: "  "4121"
## [1] "s: "  "4122"
## [1] "s: "  "4123"
## [1] "s: "  "4124"
## [1] "s: "  "4125"
## [1] "s: "  "4126"
## [1] "s: "  "4127"
## [1] "s: "  "4128"
## [1] "s: "  "4129"
## [1] "s: "  "4130"
## [1] "s: "  "4131"
## [1] "s: "  "4132"
## [1] "s: "  "4133"
## [1] "s: "  "4134"
## [1] "s: "  "4135"
## [1] "s: "  "4136"
## [1] "s: "  "4137"
## [1] "s: "  "4138"
## [1] "s: "  "4139"
## [1] "s: "  "4140"
## [1] "s: "  "4141"
## [1] "s: "  "4142"
## [1] "s: "  "4143"
## [1] "s: "  "4144"
## [1] "s: "  "4145"
## [1] "s: "  "4146"
## [1] "s: "  "4147"
## [1] "s: "  "4148"
## [1] "s: "  "4149"
## [1] "s: "  "4150"
## [1] "s: "  "4151"
## [1] "s: "  "4152"
## [1] "s: "  "4153"
## [1] "s: "  "4154"
## [1] "s: "  "4155"
## [1] "s: "  "4156"
## [1] "s: "  "4157"
## [1] "s: "  "4158"
## [1] "s: "  "4159"
## [1] "s: "  "4160"
## [1] "s: "  "4161"
## [1] "s: "  "4162"
## [1] "s: "  "4163"
## [1] "s: "  "4164"
## [1] "s: "  "4165"
## [1] "s: "  "4166"
## [1] "s: "  "4167"
## [1] "s: "  "4168"
## [1] "s: "  "4169"
## [1] "s: "  "4170"
## [1] "s: "  "4171"
## [1] "s: "  "4172"
## [1] "s: "  "4173"
## [1] "s: "  "4174"
## [1] "s: "  "4175"
## [1] "s: "  "4176"
## [1] "s: "  "4177"
## [1] "s: "  "4178"
## [1] "s: "  "4179"
## [1] "s: "  "4180"
## [1] "s: "  "4181"
## [1] "s: "  "4182"
## [1] "s: "  "4183"
## [1] "s: "  "4184"
## [1] "s: "  "4185"
## [1] "s: "  "4186"
## [1] "s: "  "4187"
## [1] "s: "  "4188"
## [1] "s: "  "4189"
## [1] "s: "  "4190"
## [1] "s: "  "4191"
## [1] "s: "  "4192"
## [1] "s: "  "4193"
## [1] "s: "  "4194"
## [1] "s: "  "4195"
## [1] "s: "  "4196"
## [1] "s: "  "4197"
## [1] "s: "  "4198"
## [1] "s: "  "4199"
## [1] "s: "  "4200"
## [1] "s: "  "4201"
## [1] "s: "  "4202"
## [1] "s: "  "4203"
## [1] "s: "  "4204"
## [1] "s: "  "4205"
## [1] "s: "  "4206"
## [1] "s: "  "4207"
## [1] "s: "  "4208"
## [1] "s: "  "4209"
## [1] "s: "  "4210"
## [1] "s: "  "4211"
## [1] "s: "  "4212"
## [1] "s: "  "4213"
## [1] "s: "  "4214"
## [1] "s: "  "4215"
## [1] "s: "  "4216"
## [1] "s: "  "4217"
## [1] "s: "  "4218"
## [1] "s: "  "4219"
## [1] "s: "  "4220"
## [1] "s: "  "4221"
## [1] "s: "  "4222"
## [1] "s: "  "4223"
## [1] "s: "  "4224"
## [1] "s: "  "4225"
## [1] "s: "  "4226"
## [1] "s: "  "4227"
## [1] "s: "  "4228"
## [1] "s: "  "4229"
## [1] "s: "  "4230"
## [1] "s: "  "4231"
## [1] "s: "  "4232"
## [1] "s: "  "4233"
## [1] "s: "  "4234"
## [1] "s: "  "4235"
## [1] "s: "  "4236"
## [1] "s: "  "4237"
## [1] "s: "  "4238"
## [1] "s: "  "4239"
## [1] "s: "  "4240"
## [1] "s: "  "4241"
## [1] "s: "  "4242"
## [1] "s: "  "4243"
## [1] "s: "  "4244"
## [1] "s: "  "4245"
## [1] "s: "  "4246"
## [1] "s: "  "4247"
## [1] "s: "  "4248"
## [1] "s: "  "4249"
## [1] "s: "  "4250"
## [1] "s: "  "4251"
## [1] "s: "  "4252"
## [1] "s: "  "4253"
## [1] "s: "  "4254"
## [1] "s: "  "4255"
## [1] "s: "  "4256"
## [1] "s: "  "4257"
## [1] "s: "  "4258"
## [1] "s: "  "4259"
## [1] "s: "  "4260"
## [1] "s: "  "4261"
## [1] "s: "  "4262"
## [1] "s: "  "4263"
## [1] "s: "  "4264"
## [1] "s: "  "4265"
## [1] "s: "  "4266"
## [1] "s: "  "4267"
## [1] "s: "  "4268"
## [1] "s: "  "4269"
## [1] "s: "  "4270"
## [1] "s: "  "4271"
## [1] "s: "  "4272"
## [1] "s: "  "4273"
## [1] "s: "  "4274"
## [1] "s: "  "4275"
## [1] "s: "  "4276"
## [1] "s: "  "4277"
## [1] "s: "  "4278"
## [1] "s: "  "4279"
## [1] "s: "  "4280"
## [1] "s: "  "4281"
## [1] "s: "  "4282"
## [1] "s: "  "4283"
## [1] "s: "  "4284"
## [1] "s: "  "4285"
## [1] "s: "  "4286"
## [1] "s: "  "4287"
## [1] "s: "  "4288"
## [1] "s: "  "4289"
## [1] "s: "  "4290"
## [1] "s: "  "4291"
## [1] "s: "  "4292"
## [1] "s: "  "4293"
## [1] "s: "  "4294"
## [1] "s: "  "4295"
## [1] "s: "  "4296"
## [1] "s: "  "4297"
## [1] "s: "  "4298"
## [1] "s: "  "4299"
## [1] "s: "  "4300"
## [1] "s: "  "4301"
## [1] "s: "  "4302"
## [1] "s: "  "4303"
## [1] "s: "  "4304"
## [1] "s: "  "4305"
## [1] "s: "  "4306"
## [1] "s: "  "4307"
## [1] "s: "  "4308"
## [1] "s: "  "4309"
## [1] "s: "  "4310"
## [1] "s: "  "4311"
## [1] "s: "  "4312"
## [1] "s: "  "4313"
## [1] "s: "  "4314"
## [1] "s: "  "4315"
## [1] "s: "  "4316"
## [1] "s: "  "4317"
## [1] "s: "  "4318"
## [1] "s: "  "4319"
## [1] "s: "  "4320"
## [1] "s: "  "4321"
## [1] "s: "  "4322"
## [1] "s: "  "4323"
## [1] "s: "  "4324"
## [1] "s: "  "4325"
## [1] "s: "  "4326"
## [1] "s: "  "4327"
## [1] "s: "  "4328"
## [1] "s: "  "4329"
## [1] "s: "  "4330"
## [1] "s: "  "4331"
## [1] "s: "  "4332"
## [1] "s: "  "4333"
## [1] "s: "  "4334"
## [1] "s: "  "4335"
## [1] "s: "  "4336"
## [1] "s: "  "4337"
## [1] "s: "  "4338"
## [1] "s: "  "4339"
## [1] "s: "  "4340"
## [1] "s: "  "4341"
## [1] "s: "  "4342"
## [1] "s: "  "4343"
## [1] "s: "  "4344"
## [1] "s: "  "4345"
## [1] "s: "  "4346"
## [1] "s: "  "4347"
## [1] "s: "  "4348"
## [1] "s: "  "4349"
## [1] "s: "  "4350"
## [1] "s: "  "4351"
## [1] "s: "  "4352"
## [1] "s: "  "4353"
## [1] "s: "  "4354"
## [1] "s: "  "4355"
## [1] "s: "  "4356"
## [1] "s: "  "4357"
## [1] "s: "  "4358"
## [1] "s: "  "4359"
## [1] "s: "  "4360"
## [1] "s: "  "4361"
## [1] "s: "  "4362"
## [1] "s: "  "4363"
## [1] "s: "  "4364"
## [1] "s: "  "4365"
## [1] "s: "  "4366"
## [1] "s: "  "4367"
## [1] "s: "  "4368"
## [1] "s: "  "4369"
## [1] "s: "  "4370"
## [1] "s: "  "4371"
## [1] "s: "  "4372"
## [1] "s: "  "4373"
## [1] "s: "  "4374"
## [1] "s: "  "4375"
## [1] "s: "  "4376"
## [1] "s: "  "4377"
## [1] "s: "  "4378"
## [1] "s: "  "4379"
## [1] "s: "  "4380"
## [1] "s: "  "4381"
## [1] "s: "  "4382"
## [1] "s: "  "4383"
## [1] "s: "  "4384"
## [1] "s: "  "4385"
## [1] "s: "  "4386"
## [1] "s: "  "4387"
## [1] "s: "  "4388"
## [1] "s: "  "4389"
## [1] "s: "  "4390"
## [1] "s: "  "4391"
## [1] "s: "  "4392"
## [1] "s: "  "4393"
## [1] "s: "  "4394"
## [1] "s: "  "4395"
## [1] "s: "  "4396"
## [1] "s: "  "4397"
## [1] "s: "  "4398"
## [1] "s: "  "4399"
## [1] "s: "  "4400"
## [1] "s: "  "4401"
## [1] "s: "  "4402"
## [1] "s: "  "4403"
## [1] "s: "  "4404"
## [1] "s: "  "4405"
## [1] "s: "  "4406"
## [1] "s: "  "4407"
## [1] "s: "  "4408"
## [1] "s: "  "4409"
## [1] "s: "  "4410"
## [1] "s: "  "4411"
## [1] "s: "  "4412"
## [1] "s: "  "4413"
## [1] "s: "  "4414"
## [1] "s: "  "4415"
## [1] "s: "  "4416"
## [1] "s: "  "4417"
## [1] "s: "  "4418"
## [1] "s: "  "4419"
## [1] "s: "  "4420"
## [1] "s: "  "4421"
## [1] "s: "  "4422"
## [1] "s: "  "4423"
## [1] "s: "  "4424"
## [1] "s: "  "4425"
## [1] "s: "  "4426"
## [1] "s: "  "4427"
## [1] "s: "  "4428"
## [1] "s: "  "4429"
## [1] "s: "  "4430"
## [1] "s: "  "4431"
## [1] "s: "  "4432"
## [1] "s: "  "4433"
## [1] "s: "  "4434"
## [1] "s: "  "4435"
## [1] "s: "  "4436"
## [1] "s: "  "4437"
## [1] "s: "  "4438"
## [1] "s: "  "4439"
## [1] "s: "  "4440"
## [1] "s: "  "4441"
## [1] "s: "  "4442"
## [1] "s: "  "4443"
## [1] "s: "  "4444"
## [1] "s: "  "4445"
## [1] "s: "  "4446"
## [1] "s: "  "4447"
## [1] "s: "  "4448"
## [1] "s: "  "4449"
## [1] "s: "  "4450"
## [1] "s: "  "4451"
## [1] "s: "  "4452"
## [1] "s: "  "4453"
## [1] "s: "  "4454"
## [1] "s: "  "4455"
## [1] "s: "  "4456"
## [1] "s: "  "4457"
## [1] "s: "  "4458"
## [1] "s: "  "4459"
## [1] "s: "  "4460"
## [1] "s: "  "4461"
## [1] "s: "  "4462"
## [1] "s: "  "4463"
## [1] "s: "  "4464"
## [1] "s: "  "4465"
## [1] "s: "  "4466"
## [1] "s: "  "4467"
## [1] "s: "  "4468"
## [1] "s: "  "4469"
## [1] "s: "  "4470"
## [1] "s: "  "4471"
## [1] "s: "  "4472"
## [1] "s: "  "4473"
## [1] "s: "  "4474"
## [1] "s: "  "4475"
## [1] "s: "  "4476"
## [1] "s: "  "4477"
## [1] "s: "  "4478"
## [1] "s: "  "4479"
## [1] "s: "  "4480"
## [1] "s: "  "4481"
## [1] "s: "  "4482"
## [1] "s: "  "4483"
## [1] "s: "  "4484"
## [1] "s: "  "4485"
## [1] "s: "  "4486"
## [1] "s: "  "4487"
## [1] "s: "  "4488"
## [1] "s: "  "4489"
## [1] "s: "  "4490"
## [1] "s: "  "4491"
## [1] "s: "  "4492"
## [1] "s: "  "4493"
## [1] "s: "  "4494"
## [1] "s: "  "4495"
## [1] "s: "  "4496"
## [1] "s: "  "4497"
## [1] "s: "  "4498"
## [1] "s: "  "4499"
## [1] "s: "  "4500"
## [1] "s: "  "4501"
## [1] "s: "  "4502"
## [1] "s: "  "4503"
## [1] "s: "  "4504"
## [1] "s: "  "4505"
## [1] "s: "  "4506"
## [1] "s: "  "4507"
## [1] "s: "  "4508"
## [1] "s: "  "4509"
## [1] "s: "  "4510"
## [1] "s: "  "4511"
## [1] "s: "  "4512"
## [1] "s: "  "4513"
## [1] "s: "  "4514"
## [1] "s: "  "4515"
## [1] "s: "  "4516"
## [1] "s: "  "4517"
## [1] "s: "  "4518"
## [1] "s: "  "4519"
## [1] "s: "  "4520"
## [1] "s: "  "4521"
## [1] "s: "  "4522"
## [1] "s: "  "4523"
## [1] "s: "  "4524"
## [1] "s: "  "4525"
## [1] "s: "  "4526"
## [1] "s: "  "4527"
## [1] "s: "  "4528"
## [1] "s: "  "4529"
## [1] "s: "  "4530"
## [1] "s: "  "4531"
## [1] "s: "  "4532"
## [1] "s: "  "4533"
## [1] "s: "  "4534"
## [1] "s: "  "4535"
## [1] "s: "  "4536"
## [1] "s: "  "4537"
## [1] "s: "  "4538"
## [1] "s: "  "4539"
## [1] "s: "  "4540"
## [1] "s: "  "4541"
## [1] "s: "  "4542"
## [1] "s: "  "4543"
## [1] "s: "  "4544"
## [1] "s: "  "4545"
## [1] "s: "  "4546"
## [1] "s: "  "4547"
## [1] "s: "  "4548"
## [1] "s: "  "4549"
## [1] "s: "  "4550"
## [1] "s: "  "4551"
## [1] "s: "  "4552"
## [1] "s: "  "4553"
## [1] "s: "  "4554"
## [1] "s: "  "4555"
## [1] "s: "  "4556"
## [1] "s: "  "4557"
## [1] "s: "  "4558"
## [1] "s: "  "4559"
## [1] "s: "  "4560"
## [1] "s: "  "4561"
## [1] "s: "  "4562"
## [1] "s: "  "4563"
## [1] "s: "  "4564"
## [1] "s: "  "4565"
## [1] "s: "  "4566"
## [1] "s: "  "4567"
## [1] "s: "  "4568"
## [1] "s: "  "4569"
## [1] "s: "  "4570"
## [1] "s: "  "4571"
## [1] "s: "  "4572"
## [1] "s: "  "4573"
## [1] "s: "  "4574"
## [1] "s: "  "4575"
## [1] "s: "  "4576"
## [1] "s: "  "4577"
## [1] "s: "  "4578"
## [1] "s: "  "4579"
## [1] "s: "  "4580"
## [1] "s: "  "4581"
## [1] "s: "  "4582"
## [1] "s: "  "4583"
## [1] "s: "  "4584"
## [1] "s: "  "4585"
## [1] "s: "  "4586"
## [1] "s: "  "4587"
## [1] "s: "  "4588"
## [1] "s: "  "4589"
## [1] "s: "  "4590"
## [1] "s: "  "4591"
## [1] "s: "  "4592"
## [1] "s: "  "4593"
## [1] "s: "  "4594"
## [1] "s: "  "4595"
## [1] "s: "  "4596"
## [1] "s: "  "4597"
## [1] "s: "  "4598"
## [1] "s: "  "4599"
## [1] "s: "  "4600"
## [1] "s: "  "4601"
## [1] "s: "  "4602"
## [1] "s: "  "4603"
## [1] "s: "  "4604"
## [1] "s: "  "4605"
## [1] "s: "  "4606"
## [1] "s: "  "4607"
## [1] "s: "  "4608"
## [1] "s: "  "4609"
## [1] "s: "  "4610"
## [1] "s: "  "4611"
## [1] "s: "  "4612"
## [1] "s: "  "4613"
## [1] "s: "  "4614"
## [1] "s: "  "4615"
## [1] "s: "  "4616"
## [1] "s: "  "4617"
## [1] "s: "  "4618"
## [1] "s: "  "4619"
## [1] "s: "  "4620"
## [1] "s: "  "4621"
## [1] "s: "  "4622"
## [1] "s: "  "4623"
## [1] "s: "  "4624"
## [1] "s: "  "4625"
## [1] "s: "  "4626"
## [1] "s: "  "4627"
## [1] "s: "  "4628"
## [1] "s: "  "4629"
## [1] "s: "  "4630"
## [1] "s: "  "4631"
## [1] "s: "  "4632"
## [1] "s: "  "4633"
## [1] "s: "  "4634"
## [1] "s: "  "4635"
## [1] "s: "  "4636"
## [1] "s: "  "4637"
## [1] "s: "  "4638"
## [1] "s: "  "4639"
## [1] "s: "  "4640"
## [1] "s: "  "4641"
## [1] "s: "  "4642"
## [1] "s: "  "4643"
## [1] "s: "  "4644"
## [1] "s: "  "4645"
## [1] "s: "  "4646"
## [1] "s: "  "4647"
## [1] "s: "  "4648"
## [1] "s: "  "4649"
## [1] "s: "  "4650"
## [1] "s: "  "4651"
## [1] "s: "  "4652"
## [1] "s: "  "4653"
## [1] "s: "  "4654"
## [1] "s: "  "4655"
## [1] "s: "  "4656"
## [1] "s: "  "4657"
## [1] "s: "  "4658"
## [1] "s: "  "4659"
## [1] "s: "  "4660"
## [1] "s: "  "4661"
## [1] "s: "  "4662"
## [1] "s: "  "4663"
## [1] "s: "  "4664"
## [1] "s: "  "4665"
## [1] "s: "  "4666"
## [1] "s: "  "4667"
## [1] "s: "  "4668"
## [1] "s: "  "4669"
## [1] "s: "  "4670"
## [1] "s: "  "4671"
## [1] "s: "  "4672"
## [1] "s: "  "4673"
## [1] "s: "  "4674"
## [1] "s: "  "4675"
## [1] "s: "  "4676"
## [1] "s: "  "4677"
## [1] "s: "  "4678"
## [1] "s: "  "4679"
## [1] "s: "  "4680"
## [1] "s: "  "4681"
## [1] "s: "  "4682"
## [1] "s: "  "4683"
## [1] "s: "  "4684"
## [1] "s: "  "4685"
## [1] "s: "  "4686"
## [1] "s: "  "4687"
## [1] "s: "  "4688"
## [1] "s: "  "4689"
## [1] "s: "  "4690"
## [1] "s: "  "4691"
## [1] "s: "  "4692"
## [1] "s: "  "4693"
## [1] "s: "  "4694"
## [1] "s: "  "4695"
## [1] "s: "  "4696"
## [1] "s: "  "4697"
## [1] "s: "  "4698"
## [1] "s: "  "4699"
## [1] "s: "  "4700"
## [1] "s: "  "4701"
## [1] "s: "  "4702"
## [1] "s: "  "4703"
## [1] "s: "  "4704"
## [1] "s: "  "4705"
## [1] "s: "  "4706"
## [1] "s: "  "4707"
## [1] "s: "  "4708"
## [1] "s: "  "4709"
## [1] "s: "  "4710"
## [1] "s: "  "4711"
## [1] "s: "  "4712"
## [1] "s: "  "4713"
## [1] "s: "  "4714"
## [1] "s: "  "4715"
## [1] "s: "  "4716"
## [1] "s: "  "4717"
## [1] "s: "  "4718"
## [1] "s: "  "4719"
## [1] "s: "  "4720"
## [1] "s: "  "4721"
## [1] "s: "  "4722"
## [1] "s: "  "4723"
## [1] "s: "  "4724"
## [1] "s: "  "4725"
## [1] "s: "  "4726"
## [1] "s: "  "4727"
## [1] "s: "  "4728"
## [1] "s: "  "4729"
## [1] "s: "  "4730"
## [1] "s: "  "4731"
## [1] "s: "  "4732"
## [1] "s: "  "4733"
## [1] "s: "  "4734"
## [1] "s: "  "4735"
## [1] "s: "  "4736"
## [1] "s: "  "4737"
## [1] "s: "  "4738"
## [1] "s: "  "4739"
## [1] "s: "  "4740"
## [1] "s: "  "4741"
## [1] "s: "  "4742"
## [1] "s: "  "4743"
## [1] "s: "  "4744"
## [1] "s: "  "4745"
## [1] "s: "  "4746"
## [1] "s: "  "4747"
## [1] "s: "  "4748"
## [1] "s: "  "4749"
## [1] "s: "  "4750"
## [1] "s: "  "4751"
## [1] "s: "  "4752"
## [1] "s: "  "4753"
## [1] "s: "  "4754"
## [1] "s: "  "4755"
## [1] "s: "  "4756"
## [1] "s: "  "4757"
## [1] "s: "  "4758"
## [1] "s: "  "4759"
## [1] "s: "  "4760"
## [1] "s: "  "4761"
## [1] "s: "  "4762"
## [1] "s: "  "4763"
## [1] "s: "  "4764"
## [1] "s: "  "4765"
## [1] "s: "  "4766"
## [1] "s: "  "4767"
## [1] "s: "  "4768"
## [1] "s: "  "4769"
## [1] "s: "  "4770"
## [1] "s: "  "4771"
## [1] "s: "  "4772"
## [1] "s: "  "4773"
## [1] "s: "  "4774"
## [1] "s: "  "4775"
## [1] "s: "  "4776"
## [1] "s: "  "4777"
## [1] "s: "  "4778"
## [1] "s: "  "4779"
## [1] "s: "  "4780"
## [1] "s: "  "4781"
## [1] "s: "  "4782"
## [1] "s: "  "4783"
## [1] "s: "  "4784"
## [1] "s: "  "4785"
## [1] "s: "  "4786"
## [1] "s: "  "4787"
## [1] "s: "  "4788"
## [1] "s: "  "4789"
## [1] "s: "  "4790"
## [1] "s: "  "4791"
## [1] "s: "  "4792"
## [1] "s: "  "4793"
## [1] "s: "  "4794"
## [1] "s: "  "4795"
## [1] "s: "  "4796"
## [1] "s: "  "4797"
## [1] "s: "  "4798"
## [1] "s: "  "4799"
## [1] "s: "  "4800"
## [1] "s: "  "4801"
## [1] "s: "  "4802"
## [1] "s: "  "4803"
## [1] "s: "  "4804"
## [1] "s: "  "4805"
## [1] "s: "  "4806"
## [1] "s: "  "4807"
## [1] "s: "  "4808"
## [1] "s: "  "4809"
## [1] "s: "  "4810"
## [1] "s: "  "4811"
## [1] "s: "  "4812"
## [1] "s: "  "4813"
## [1] "s: "  "4814"
## [1] "s: "  "4815"
## [1] "s: "  "4816"
## [1] "s: "  "4817"
## [1] "s: "  "4818"
## [1] "s: "  "4819"
## [1] "s: "  "4820"
## [1] "s: "  "4821"
## [1] "s: "  "4822"
## [1] "s: "  "4823"
## [1] "s: "  "4824"
## [1] "s: "  "4825"
## [1] "s: "  "4826"
## [1] "s: "  "4827"
## [1] "s: "  "4828"
## [1] "s: "  "4829"
## [1] "s: "  "4830"
## [1] "s: "  "4831"
## [1] "s: "  "4832"
## [1] "s: "  "4833"
## [1] "s: "  "4834"
## [1] "s: "  "4835"
## [1] "s: "  "4836"
## [1] "s: "  "4837"
## [1] "s: "  "4838"
## [1] "s: "  "4839"
## [1] "s: "  "4840"
## [1] "s: "  "4841"
## [1] "s: "  "4842"
## [1] "s: "  "4843"
## [1] "s: "  "4844"
## [1] "s: "  "4845"
## [1] "s: "  "4846"
## [1] "s: "  "4847"
## [1] "s: "  "4848"
## [1] "s: "  "4849"
## [1] "s: "  "4850"
## [1] "s: "  "4851"
## [1] "s: "  "4852"
## [1] "s: "  "4853"
## [1] "s: "  "4854"
## [1] "s: "  "4855"
## [1] "s: "  "4856"
## [1] "s: "  "4857"
## [1] "s: "  "4858"
## [1] "s: "  "4859"
## [1] "s: "  "4860"
## [1] "s: "  "4861"
## [1] "s: "  "4862"
## [1] "s: "  "4863"
## [1] "s: "  "4864"
## [1] "s: "  "4865"
## [1] "s: "  "4866"
## [1] "s: "  "4867"
## [1] "s: "  "4868"
## [1] "s: "  "4869"
## [1] "s: "  "4870"
## [1] "s: "  "4871"
## [1] "s: "  "4872"
## [1] "s: "  "4873"
## [1] "s: "  "4874"
## [1] "s: "  "4875"
## [1] "s: "  "4876"
## [1] "s: "  "4877"
## [1] "s: "  "4878"
## [1] "s: "  "4879"
## [1] "s: "  "4880"
## [1] "s: "  "4881"
## [1] "s: "  "4882"
## [1] "s: "  "4883"
## [1] "s: "  "4884"
## [1] "s: "  "4885"
## [1] "s: "  "4886"
## [1] "s: "  "4887"
## [1] "s: "  "4888"
## [1] "s: "  "4889"
## [1] "s: "  "4890"
## [1] "s: "  "4891"
## [1] "s: "  "4892"
## [1] "s: "  "4893"
## [1] "s: "  "4894"
## [1] "s: "  "4895"
## [1] "s: "  "4896"
## [1] "s: "  "4897"
## [1] "s: "  "4898"
## [1] "s: "  "4899"
## [1] "s: "  "4900"
## [1] "s: "  "4901"
## [1] "s: "  "4902"
## [1] "s: "  "4903"
## [1] "s: "  "4904"
## [1] "s: "  "4905"
## [1] "s: "  "4906"
## [1] "s: "  "4907"
## [1] "s: "  "4908"
## [1] "s: "  "4909"
## [1] "s: "  "4910"
## [1] "s: "  "4911"
## [1] "s: "  "4912"
## [1] "s: "  "4913"
## [1] "s: "  "4914"
## [1] "s: "  "4915"
## [1] "s: "  "4916"
## [1] "s: "  "4917"
## [1] "s: "  "4918"
## [1] "s: "  "4919"
## [1] "s: "  "4920"
## [1] "s: "  "4921"
## [1] "s: "  "4922"
## [1] "s: "  "4923"
## [1] "s: "  "4924"
## [1] "s: "  "4925"
## [1] "s: "  "4926"
## [1] "s: "  "4927"
## [1] "s: "  "4928"
## [1] "s: "  "4929"
## [1] "s: "  "4930"
## [1] "s: "  "4931"
## [1] "s: "  "4932"
## [1] "s: "  "4933"
## [1] "s: "  "4934"
## [1] "s: "  "4935"
## [1] "s: "  "4936"
## [1] "s: "  "4937"
## [1] "s: "  "4938"
## [1] "s: "  "4939"
## [1] "s: "  "4940"
## [1] "s: "  "4941"
## [1] "s: "  "4942"
## [1] "s: "  "4943"
## [1] "s: "  "4944"
## [1] "s: "  "4945"
## [1] "s: "  "4946"
## [1] "s: "  "4947"
## [1] "s: "  "4948"
## [1] "s: "  "4949"
## [1] "s: "  "4950"
## [1] "s: "  "4951"
## [1] "s: "  "4952"
## [1] "s: "  "4953"
## [1] "s: "  "4954"
## [1] "s: "  "4955"
## [1] "s: "  "4956"
## [1] "s: "  "4957"
## [1] "s: "  "4958"
## [1] "s: "  "4959"
## [1] "s: "  "4960"
## [1] "s: "  "4961"
## [1] "s: "  "4962"
## [1] "s: "  "4963"
## [1] "s: "  "4964"
## [1] "s: "  "4965"
## [1] "s: "  "4966"
## [1] "s: "  "4967"
## [1] "s: "  "4968"
## [1] "s: "  "4969"
## [1] "s: "  "4970"
## [1] "s: "  "4971"
## [1] "s: "  "4972"
## [1] "s: "  "4973"
## [1] "s: "  "4974"
## [1] "s: "  "4975"
## [1] "s: "  "4976"
## [1] "s: "  "4977"
## [1] "s: "  "4978"
## [1] "s: "  "4979"
## [1] "s: "  "4980"
## [1] "s: "  "4981"
## [1] "s: "  "4982"
## [1] "s: "  "4983"
## [1] "s: "  "4984"
## [1] "s: "  "4985"
## [1] "s: "  "4986"
## [1] "s: "  "4987"
## [1] "s: "  "4988"
## [1] "s: "  "4989"
## [1] "s: "  "4990"
## [1] "s: "  "4991"
## [1] "s: "  "4992"
## [1] "s: "  "4993"
## [1] "s: "  "4994"
## [1] "s: "  "4995"
## [1] "s: "  "4996"
## [1] "s: "  "4997"
## [1] "s: "  "4998"
## [1] "s: "  "4999"
## [1] "s: "  "5000"
pred_mean_healthy <- unlist(map(MEAN, mean))

pred_df_healthy = tibble(
  time = t,
  pred_rate = pred_mean_healthy,
  emp_rate = y
)
pred_df_healthy %>% 
  ggplot() +
  geom_line(aes(x = time, y = emp_rate), col = "blue") +
  geom_line(aes(x = time, y = pred_rate), col = "orange")

par(mfrow = c(3, 2))
for(j in sample(12, 3)){
  plot(ALPHA[[j]], type = 'l', ylab = paste("Alpha", j), xlab = "Iter", col = "purple",
       main = paste("Acceptance ratio =", round(length(unique(ALPHA[[j]])) / S, 3)))
  acf(ALPHA[[j]], ylab = paste("Alpha", j), main = paste("Series Alpha", j))
}

par(mfrow = c(3, 2))
for(j in sample(17, 3)){
  plot(DELTA[[j]], type = 'l', ylab = paste("Delta", j), xlab = "Iter", col = "purple",
       main = paste("Acceptance ratio =", round(length(unique(DELTA[[j]])) / S, 3)))
  acf(DELTA[[j]], ylab = paste("Delta", j), main = paste("Series Delta", j))
}

par(mfrow = c(1,2))
plot(DC, type = 'l', ylab = "DC", xlab = "Iter", col = "purple",
       main = paste("Acceptance ratio =", round(length(unique(DC)) / S, 3)))
acf(DC, ylab = "DC", main = "Series DC")

par(mfrow = c(1,2))
plot(TAO, type = 'l', ylab = "Tao", xlab = "Iter", col = "purple",
    main = paste("Acceptance ratio =", round(length(unique(TAO)) / S, 3)))
acf(TAO, ylab = "Tao", main = "Series Tao")

alpha_names <- map_chr(1:12, function(x) str_c("alpha_", x, sep = ""))
alpha_df <- map_dfc(ALPHA, function(x) x)
colnames(alpha_df) <- alpha_names

alpha_df <- alpha_df %>%
  pivot_longer(1:12, names_to = "index", values_to = "alpha_1") %>%
  mutate(index = factor(index, levels = alpha_names))

colnames(alpha_df)[2] <- "value"

alpha_df %>%
  filter(index %in% alpha_names[1:6]) %>%
  ggplot() +
  geom_density(aes(x = value, fill = index), alpha = 0.7) 

alpha_df %>%
  filter(index %in% alpha_names[7:12]) %>%
  ggplot() +
  geom_density(aes(x = value, fill = index), alpha = 0.7) 

delta_names <- map_chr(1:17, function(x) str_c("delta_", x, sep = ""))
delta_df <- map_dfc(DELTA, function(x) x)
colnames(delta_df) <- delta_names

delta_df <- delta_df %>%
  pivot_longer(1:17, names_to = "index", values_to = "delta_1") %>%
  mutate(index = factor(index, levels = delta_names))

colnames(delta_df)[2] <- "value"

delta_df %>%
  filter(index %in% delta_names[1:6]) %>%
  ggplot() +
  geom_density(aes(x = value, fill = index), alpha = 0.7) 

delta_df %>%
  filter(index %in% delta_names[7:12]) %>%
  ggplot() +
  geom_density(aes(x = value, fill = index), alpha = 0.7) 

delta_df %>%
  filter(index %in% delta_names[13:17]) %>%
  ggplot() +
  geom_density(aes(x = value, fill = index), alpha = 0.7) 

for(j in sample(length(MEAN), 6)){
  plot_2d <- ggplot() +
    geom_point(aes(x = MEAN[[j]], y = TAO)) +
    geom_density_2d(aes(x = MEAN[[j]], y = TAO))
  
  print(plot_2d)
}

samples_arr <- read_csv("samples.csv", skip = 1)
## Parsed with column specification:
## cols(
##   `'seconds'` = col_double(),
##   `'mV'` = col_double()
## )
colnames(samples_arr) <- c("time", "signal")
y <- samples_arr$signal[76:662] #signal 0.208 max(samples$signal[0:250]) which(samples$signal[0:250] == 0.58)
t <- samples_arr$time[76:662] #signal 0.705 max(samples$signal[0:700]) which(samples$signal[0:700] == 0.705)
S <- 5000
burn <- 5000
beta_curr <- 1
alpha_curr <- rep(1, 12)
delta_curr<- seq(0, max(t), length.out = 17) 
end_time <- max(t)
tao_curr <- 1
mean_curr <- rep(mean(y), length(y)) #mean_curr is a vector cause mean_star of y is a vector so I think r only lets me run it if it is a vector
dc_curr <- 0.5
miu_curr <- mean_curr - dc_curr

#storage

DC <- NULL
ALPHA <- rep(list(NULL), 12)
BETA <- NULL
DELTA <- rep(list(NULL), 17)
TAO <- NULL
MEAN <- rep(list(NULL), length(y))
for(s in 1:S) {
  #print(c("s: ", s))
  
  #1 draw of Beta
  
  beta_lst <- beta_one_samp(beta_curr, BETA, s)
  beta_curr <- beta_lst[[1]]
  BETA <- beta_lst[[2]]
  #print(c("beta: ", beta_curr))
  
  #1 draw of Alphas (1-12)
  
  for (k in 1:12) {
    alpha_lst <- alpha_one_samp(k, alpha_curr[k], ALPHA[[k]], s)
    alpha_curr[k] <- alpha_lst[[1]]
    ALPHA[[k]] <- alpha_lst[[2]]
  }
  
  #print(c("ALPHA: ", alpha_curr))
  
  #1 draw of Deltas (1-17)
  
  delta_lst <- delta_one_samp(delta_curr[1], 0, delta_curr[2], DELTA[[1]], s)
  delta_curr[1] <- delta_lst[[1]]
  DELTA[[1]] <- delta_lst[[2]]
  
  for (k in 2:16) {
    delta_lst <- delta_one_samp(delta_curr[k], delta_curr[k-1], delta_curr[k+1], DELTA[[k]], s)
    delta_curr[k] <- delta_lst[[1]]
    DELTA[[k]] <- delta_lst[[2]]
  }
  
  delta_lst <- delta_one_samp(delta_curr[17], delta_curr[[16]], end_time,  DELTA[[17]], s)
  delta_curr[17] <- delta_lst[[1]]
  DELTA[[17]] <- delta_lst[[2]]
  
  #print(c("DELTA: ", delta_curr))
  
  
  #1 draw of Tao 
  tao_lst <- tao_one_samp_cond(tao_curr, TAO, s)
  tao_curr <- tao_lst[[1]]
  TAO <- tao_lst[[2]]

  #print(c("TAO: ", tao_curr))
  
  
  #1 draw of DC 
  dc_lst <- dc_one_samp(dc_curr, DC, s)
  dc_curr <- dc_lst[[1]]
  DC <- dc_lst[[2]]
  
  #print(c("DC: ", dc_curr))
  
  mean_curr <- dc_curr + gen_miu_star(t, alpha_curr = alpha_curr, beta_curr = beta_curr, delta_curr = delta_curr)
  for (k in 1:length(y)) {
    MEAN[[k]] <- c(MEAN[[k]], mean_curr[[k]])
  }
  
  #print(c("MEAN CURR: ", mean_curr))
}
pred_mean_arr <- unlist(map(MEAN, mean))

pred_df_arr = tibble(
  time = t,
  pred_rate = pred_mean_arr,
  emp_rate = y
)
pred_df_arr %>% 
  ggplot() +
  geom_line(aes(x = time, y = emp_rate), col = "blue") +
  geom_line(aes(x = time, y = pred_rate), col = "orange")

par(mfrow = c(3, 2))
for(j in sample(12, 3)){
  plot(ALPHA[[j]], type = 'l', ylab = paste("Alpha", j), xlab = "Iter", col = "purple",
       main = paste("Acceptance ratio =", round(length(unique(ALPHA[[j]])) / S, 3)))
  acf(ALPHA[[j]], ylab = paste("Alpha", j), main = paste("Series Alpha", j))
}

par(mfrow = c(3, 2))
for(j in sample(17, 3)){
  plot(DELTA[[j]], type = 'l', ylab = paste("Delta", j), xlab = "Iter", col = "purple",
       main = paste("Acceptance ratio =", round(length(unique(DELTA[[j]])) / S, 3)))
  acf(DELTA[[j]], ylab = paste("Delta", j), main = paste("Series Delta", j))
}

par(mfrow = c(1,2))
plot(DC, type = 'l', ylab = "DC", xlab = "Iter", col = "purple",
       main = paste("Acceptance ratio =", round(length(unique(DC)) / S, 3)))
acf(DC, ylab = "DC", main = "Series DC")

par(mfrow = c(1,2))
plot(TAO, type = 'l', ylab = "Tao", xlab = "Iter", col = "purple",
    main = paste("Acceptance ratio =", round(length(unique(TAO)) / S, 3)))
acf(TAO, ylab = "Tao", main = "Series Tao")

alpha_names <- map_chr(1:12, function(x) str_c("alpha_", x, sep = ""))
alpha_df <- map_dfc(ALPHA, function(x) x)
colnames(alpha_df) <- alpha_names

alpha_df <- alpha_df %>%
  pivot_longer(1:12, names_to = "index", values_to = "alpha_1") %>%
  mutate(index = factor(index, levels = alpha_names))

colnames(alpha_df)[2] <- "value"

alpha_df %>%
  filter(index %in% alpha_names[1:6]) %>%
  ggplot() +
  geom_density(aes(x = value, fill = index), alpha = 0.7) 

alpha_df %>%
  filter(index %in% alpha_names[7:12]) %>%
  ggplot() +
  geom_density(aes(x = value, fill = index), alpha = 0.7) 

delta_names <- map_chr(1:17, function(x) str_c("delta_", x, sep = ""))
delta_df <- map_dfc(DELTA, function(x) x)
colnames(delta_df) <- delta_names

delta_df <- delta_df %>%
  pivot_longer(1:17, names_to = "index", values_to = "delta_1") %>%
  mutate(index = factor(index, levels = delta_names))

colnames(delta_df)[2] <- "value"

delta_df %>%
  filter(index %in% delta_names[1:6]) %>%
  ggplot() +
  geom_density(aes(x = value, fill = index), alpha = 0.7) 

delta_df %>%
  filter(index %in% delta_names[7:12]) %>%
  ggplot() +
  geom_density(aes(x = value, fill = index), alpha = 0.7) 

delta_df %>%
  filter(index %in% delta_names[13:17]) %>%
  ggplot() +
  geom_density(aes(x = value, fill = index), alpha = 0.7) 

for(j in sample(length(MEAN), 6)){
  plot_2d <- ggplot() +
    geom_point(aes(x = MEAN[[j]], y = TAO)) +
    geom_density_2d(aes(x = MEAN[[j]], y = TAO))
  
  print(plot_2d)
}